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FOREWORD 


The analytical report contained herein is submitted to the NASA Langley 
Research Center by the AiResearch Manufacturing Company of Los Angeles, 
California, in accordance with the guidelines of Paragraph 5. 7. 3. 2. I, NASA 
Statement of Work L-4947-B (Revised). 

This report describes the chemical kinetic study for the Supersonic 
Combustor Model which was performed by AiResearch in support of the combustor 
program phase of the Hypersonic Research Engine Project. 



AIRESEARCH MANUEACTURING COMPANY 
Lot Antettt, CiMorma 


UNCLASSIHED 


70-6319 
Page i i 





Sect Ion 

UNCLASSIFIED 

CONTENTS 

Page 

1 

SUMMARY 

1 -1 

2 

INTRODUCTION 

2-1 

3 

PARAMETRIC STUDIES 

3-1 


3.1 The Detailed Hydrogen-Air Chemical Kinetic Model 

3-1 


3.2 Results and Discussion 

3-2 

4 

CHEMICAL KINETICS DURING THE IGNITION DELAY PERIOD 

4-1 

5 

TWO-DIMENSIONAL EFFECTS IN A DIVERGING COMBUSTOR 

5-1 

6 

COMPARISON WITH TEST DATA 

6-1 


6.1 Results 

6-3 


6.2 Chemical Kinetic Simulation of Two-Dimensional 

6-1 1 

7 

AIM ANALYSIS 

7-1 

8 

DISCUSSION AND RECOMMENDATIONS 

8-1 


8.1 D i scuss ion 

8-1 


8.2 Future Work Recommendations 

8-1 

9 

CONCLUSIONS 

9-1 

References 

R-l 

Appendix 

A Run 238 


Appendix 

B Run 242 



70-6319 
Page i i i 


AIRESEARCM MANUFACTURING COMPANY 

Lot Anitm, 


UNCLASSIFIED 




UNCLASSIRED 


Fiqu re 

3.2- 1 

3.2- 2 

3.2- 3 

3.2- 4 

3.2- 5 

3.2- 6 

3.2- 7 

3.2- 8 

3.2- 9 

3.2- 10 

3.2- 1 I 

3.2- 12 

3.2- 13 

3.2- 14 

3.2- 15 

3.2- 16 

3.2- 17 

3.2- 18 


ILLUSTRATIONS 


of Initial Static Temperature on Combustor Length 
Effect of Initial Static Temperature on Static Pressure 

of Initial Static Temperature on Chemical Efficiency 3-5 


Species Concentration as a Function of Combustor Length 3-6 
for an Initial Static Temperature of 2200®R 

Species Concentration as a Function of Combustor Length 3-7 
for an Initial Static Temperature 

Species Concentration as a Function of Combustor Length 3-8 
for a Static Temperature of I900°R 

Effect of Initial Static Temperature on Total Ignition 3-9 

Delay 


Effect of Initial Static Pressure on Total Ignition Delay 3-13 


Initial Static Pressure on Combustor Length 3-14 

and Temperature 

Effect of Initial Static Pressure on Pressure as a 3-15 

Function of Combustor Length 

of Initial Static Pressure on Chemical Efficiency 3-16 
as a Function of Combustor Length 

Species Concentration as a Function of Combustor Length 3-17 

Effect of Initial Static Pressure on Ignition Delay 3-18 

Effect of Initial Static Pressure on Temperature as a 3-19 

Function of Combustor Length 

Effect of Initial Static Pressure on Pressure as a Func- 3-20 
tion of Combustor Length 

Effect of Initial Static Pressure on Chemical Efficiency 3-21 
as a Function of Combustor Length 

of Initial Static Pressure and Temperature on 3-22 

Ignition Delay 

Effect of Equivalence Ratio on Ignition Delay 3-24 


Page 

3-3 

3-4 


AlftESEARCM manufacturing COMPANY 

Lot AnfKM. CaMonM 


UNCLASSIFIED 


70-6319 
Page iv 




UNCLASSinED 


ILLUSTRATIONS (Continued) 


liailiie Page 

3.2- 19 Effect of Equivalence Ratio on Static Temperature 3-25 

as a Function of Combustor Length 

3.2- 20 Effect of Equivalence Ratio on Static Pressure as 3-26 

a Function of Combustor Length 

3.2- 21 Effect of Equivalence Ratio on Chemical Efficiency 3-27 

as a Function of Combustor Length 

3.2- 22 Effect of Vitiation on Ignition Delay Time 3-28 

3.2- 23 Effect of 10-Percent Carriers on Static Temperature 3-29 


as a Function of Combustor Length for Several Initial 
Static Temperatures 

3.2- 24 Effect of 10-Percent Carrier on Chemical Efficiency as a 3-30 

Function of Combustor Length for Several Initial Static 
Temperatures 

3.2- 25 Effect of 10-Percent Carriers on Static Pressure as a 3-31 

Function of Combustor Length for Several Initial 
Static Temperatures 

3.2- 26 Effect of 2-Percent Carriers on Static Temperature as a 3-32 

Function of Combustor Length for Several Initial 
Static Temperatures 

3.2- 27 Effect of 2-Percent Carriers on Chemical Efficiency as a 3-33 

Function of Combustor Length for Several Initial 
Static Temperatures 


3.2- 28 Effect of 2-Percent Carriers on Static Pressure as a 3-34 

Function of Combustor Length for Several Initial 
Static Temperatures 

3.2- 29 Effect of Percentage of Carriers on Initial Static Tern- 3-35 

peratures as a Function of Required Combustor Length and 
Chemical Efficiency 

3.2- 30 Effect of Divergence on Static Temperature as a Function 3-37 

of Combustor Length 

3.2- 31 Effect of Divergence of Chemical Efficiency as a Func- 3-38 

tion of Combustor Length 

3.2- 32 Effect of Divergence on Static Pressure as a Function 3-39 

of Combustor Length 



AIRESEARCH 


MANUFACTURING COMPANY 
Lot Anfttts, CoMormt 


UNCLASSIFIED 


70-6319 
Page v 




UNCLASSIHED 


ILLUSTRATIONS (Continued) 

Figure Page 

3.2- 33 Effect of Divergence on Static Temperature as a Func- 3-41 

tion of Combustor Length 

3.2- 34 Effect of Divergence on Chemical Efficiency as a 3-42 

Function of Combustor Length 

3.2- 35 Effect of Divergence on Static Pressure as a Function 3-43 

of Combustor Length 

3.2- 36 Effect of Divergence on Ignition Delay as a Function 3-44 

of Initial Static Temperature 

4.2- 1 Comparison of Theoretical and Experimental Ignition 4-5 

Delay Times 

5.0- 1 Mainstream Compression by Burning in the Boundary Layer 5-2 

6.1- 1 Chemical Kinetics Simulation of Cookson* s Experiment 6-4 

6.1- 2 Chemical Kinetics Simulation of Cookson* s Experiment 6-5 

6.1- 3 Chemical Kinetic Simulation of Cookson’ s Experiment 6-6 

6.1- 4 Chemical Kinetic Simulation of Cookson* s Experiment 6-7 

6.1- 5 Pressure Profiles for Cookson Combustor 6-10 

6.2- 1 Combustor Configuration Schematic 6-12 

6.2- 2 Second Stage Static Pressure Vibration in Two- 6-13 

Dimensional Combustor at T^^ = 3000®R 

6.2- 3 Second Stage Temperature Variation in Two-Dimensional 6-14 

Combustor at T^q = 3000®R 

6.2- 4 Species Concentration, Two-Dimensional Combustor at 6-15 

Tjq = 3000® R 

6.2- 5 Comparison of Calculated and Measured Static Pressures 6-16 

for Two-Dimensional Combustor Test at T-j-q = 3800®R 

6.2- 6 Static Temperature Comparison for Two-Dimensional 6-17 

Combustor at T^^ = 3800® R 

6.2- 7 Chemical Efficiency Variation at Tjq = 3800®R 6-18 

70-6319 
Page vi 


AIRESEARCH MANUFACTURING COMPANY 

Los Anfttis. CaMermt 


UNCLASSIFIED 




UNCLASSIFIED 


ILLUSTRATIONS (Continued) 


Figure Page 

6.2- 8 Species Concentration, Two-Dimensional Combustor Study 6-19 

at T^q = 3800° R 

6.2- 9 Species Concentration, Two-Dimensional Combustor Study 6-20 

at T^q = 3800°R 

6.2- 10 Mixing Models 6-22 

6.2- 11 Measured Static Pressures - Tests 238 and 242 6-23 

6.2- 12 Geometric and Calculated Areas vs Axial Position 6-24 

6.2- 13 Mixing Efficiency vs Axial Position - Run 238 6-25 

6.2- 14 Chemical Efficiency vs Axial Position 6-26 

6.2- 15 Mixing Efficiency vs Axial Position - Run 242 6-28 

6.2- 16 Geometric and Calculated Areas vs Axial Position - 6-29 

Run 242 

6.2- 17 Chemical Efficiency vs Axial Position 6-30 

7.1- 1 Static Variation at Mach 7 - AIM Combustor 7-3 

7.1- 2 Static Temperature Variation at Mach 7 - AIM Combustor 7-4 

7.1- 3 Species Concentration at Mach 7 - AIM Combustor 7-5 

7.1- 4 Species Concentration at Mach 7 - AIM Combustor 7-5 

7.1- 5 Second Stage Chemical Efficiency at Mach 7 - AIM Combustor 7-7 

7.1- 6 Static Pressure Variation at Mach 6 - AIM Combustor 7-9 

7.1- 7 Static Temperature Variation at Mach 6 - AIM Combustor 7-10 

7.1- 8 Second Stage Chemical Efficiency at Mach 6 - AIM Combustor 7-11 

7.1- 9 Species Concentration at Mach 6 - AIM Combustor 7-12 

7.1- 10 Species Concentration at Mach 6 - AIM Combustor 7-13 

7.1- 1 I Species Concentration at Mach 6 - AIM Combustor 7-14 

7.1- 12 Comparison of Species Concentration for Different Mixed 7-15 

Temperatures at Mach 6 - AIM Combustor 



AIRESEARCH MANUFACTURING COMPANY 

Lot Angdtt, CoMornta 


UNCLASSIFIED 


70-6319 
Page vi i 





UNCLASSIRED 


ILLUSTRATIONS (Continued) 


Figure Page 


7.1-13 

Temperature vs Time - Second Stage 


7-16 

7. 1-14 

Static Temperature Variation at Mach 6 - AIM Combustor 

7-18 

7. 1-15 

Effect of First Stage Efficiency on Second Stage Temperature 

7-20 

7.1-16 

Combustor Length vs Amount of First Stage Fuel 

Reaction 

7-22 

7. 1-17 

Combustor Length vs Amount of Fi rst Stage Fuel 

Reaction 

7-23 

7. 2. 1 

Static Temperature in a Constant Area Section 


7-25 

7.2 2 

Computed Pressure Profile 


7-26 

7.2-3 

Concentration of OH 


7-27 

7.2-4 

Effect of Ignition on Static Temperature 


7-28 

7.2-5 

Temperature Profile 


7-29 

7.2-6 

Gas Stream Area Expansion 


7-30 

8.2-1 

General Electric Test Data 


8-3 



AIRESEARCM MANUFAauRING COMPANY 

LotAn|t4M. CaMofma 


UNCLASSIFIED 


70-6319 
Page vi i i 



UNCLASSIFIED 


TABLES 


Table 



Page 

3.1-1 

Hydrogen-Air Reactions 


3-1 

3.1-2 

Reaction-Rate Constants 


3-10 

3.2-1 

Mass Fractions of Reacting Flow 


3-1 1 

3.2-2 

Mass Fractions of Reacting Flow with Initial 

Carriers 

3-36 

6.0-1 

Reactions Used in the Shear Layer Combustion 

Program 

6-2 

6.1-1 

Initial Conditions Assumed for Cookson' s Combustor 

6-8 



AIRESEARCH MANUFACTURING COMPANY 
Los AngMs, CsiiteKma 


UNCLASSIFIED 


70-6319 
Page ix 




UNCLASSIFIED 


I. SUMMARY 


The problems of supersonic combustion in a diverging duct were examined 
analytically, and the results of this study were compared with test results 
from the HRE two-dimensional combustor rig. The effects of temperature, pres- 
sure, duct divergence, amount of vitiation and initial concentration of 
radicals on the chemical kinetics of the hydrogen-air system were investigated. 
It was found that supersonic combustion of hydrogen-air in a diverging duct at 
low mixed temperatures (around 2000®R) is very sensitive to chemical kinetics 
effects. The ignition delay time decreased with increasing pressure but this 
effect was reversed in the pressure range of I to 15 atm depending on the 
initial temperature. The ignition process was enchanced by the presence of 
free radicals, but the relative effects diminished rapidly as the initial con- 
centration was increased. 

A detailed chemi ca 1 -kinet i c simulation of the HRE two-dimensional com- 
bustor test results indicated that low chemical efficiencies in the second 
stage were not caused by chemi ca 1 -ki net ics but were due to. inadequate mixing. 
This conclusion was further strengthed by tests reported by General Electric 
(Reference l-l) in which low chemical efficiencies were observed when fuel was 
injected into a constant-area duct followed by a divergence section, with an 
inlet static temperature as high as 2500®R. 

It is postulated that intimate mixing requires turbulence with small but 
high energy level eddies. The various injector geometries used in the two- 
dimensional combustor tests probably did not cause sufficient turbulence in 
the mainstream, evidenced by lack of sufficient flow separation. The low level 
of turbulence could be due to the favorable pressure gradient in the diverging 
duct and to the relatively high temperatures upstream of the second stage 
injectors. Further exploration in this area is urgently needed. 

The analytical simulation of the AIM engine indicated that the chemical- 
kinetics effects will not adversely affect the second-stage combustion process 
in the flight Mach number range of 6 to 8. 

The results from the two-dimensional combustor tests lead to the prediction 
of low second-stage combustor efficiencies in the AIM engine unless methods 
for enhancing the mixing process can be devised. Turbulence could be generated 
by a metallic coating on the combustor walls. This coating can be applied with 
varying degrees of surface roughness to produce different levels of turbulence. 

Currently, five additional igniter stations are being installed in the 
AIM combustor. The igniter jets are directed upstream to produce additional 
sources of turbulence as well as ignition energy. 
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2. INTRODUCTION 

The problems of supersonic combustion in a diverging duct were evident 
in the HRE two-dimensional combustor component test program. In these tests 
low chemical efficiencies, in the order of 50 to 60 percent, were encountered 
regardless of injector geometry, equivalence ratios, and inlet total tempera- 
ture. Various injector geometries and equivalence ratios aimed at improving 
the mixing, and high inlet total temperatures (3800®R) which should enhance 
the chemical reactions were tried. However, all attempts failed to improve 
the combustor performance. A subsequent survey of the literature Indicated 
that low chemical efficiencies were also obtained by a number of investigators 
(Reference |-| 1-2, and 2-1 ) whenever a diverging duct was used. Cookson 

(Reference 2-1) found that a hydrogen diffusion flame could not be maintained 
in a simple conical combustor of 0,75 deg divergence and 30 in. long. Com- 
bustion was observed to occur after the Mach disc which was formed at the 
combustor exit. Cookson's experiments strongly suggest that chemical kinetics 
have a significant affect on the problem. Therefore, an Investigation was 
initiated to study the various effects of chemical kinetics using existing 
computer programs in order to obtain a better understanding of the problem. 

The investigation was divided into four phases. First, a parametric 
study was conducted to obtain the basic understanding of the chemical kinetics 
from the effects of various independent variables such as temperature, pres- 
sure, divergence, degree of vitiation, and initial concentration of radicals; 
second, the mechanism of ignition delay was examined in a simplified analytical 
form in order to obtain the essential correlation parameters which otherwise 
could not be obtained from computer programs; and third, a comparison was made 
between the computer program results and Cookson's experiments in order to 
substantiate the validity of the kinetics programs. Fully-mixed and partially- 
mixed models were simulated in a shear-layer combustion program (Reference 2-2) 
to try and analyze and compare with the results from the two-dimensional com- 
bustor test program. Finally, a chemical kinetic analysis of the AIM engine 
was performed. A brief examination of chemical kinetics using a two-dimensional 
kinetics program was also made. 
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3. PARAMETRIC STUDIES 

An investigation of the chemical kinetics of a premixed hydrogen-air 
system was conducted. The investigation was directed toward clarifying the 
effects of static temperature, static pressure, flow velocity, equivalence 
ratio, vitiation, initial radical concentration, and combustor divergence on 
the chemical kinetics. These effects were determined through use of a one- 
dimensional chemical kinetics computer program (Reference 3-1 ). 

3.1 THE DETAILED HYDROGEN-AIR CHEMICAL KINETIC MODEL 

The reaction mechanism of the hydrogen-air reaction assumed to be com- 
prised of the II reactions is shown in Table 3.1-1. 

TABLE 3.1-1 

HYDROGEN-AIR REACTIONS 


(4) 

H2+M=H+H+M 


(9) 

02+M = 0 + 0 + M 

(chain-initiating steps) 

(II) 

H 2 + 0^ = OH + OH 


(1) 

H, + OH = H-0 + H 


(8) 

i. d 

(chain-carrying steps) 

0 + H+ M = 0H+M 

(2) 

H + O 2 = OH + 0 


(3) 

0 + H, = OH + H 


(7) 

c 

(chain-branching steps) 

H 2 O + 0 = OH + OH 

(6) 

HO 2 + «2 = H 2 O 2 + H 


(5) 

H + 0- + M = HO- + M 


(10) 

c. 2 

(chain-breaking steps) 

H + OH + M = H 2 O + M 
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The symbol M represents a "third body", i.e., any chemical species. The 
third body acts as an energy source or sink in chemical reactions. Steps 2, 
3, and 7 are branching reactions since two chain carriers appear in the 
products while only one enters into the reaction. Steps 5 and 10 serve to 
break the hydrogen-oxygen chain by removing (h) and (OH) from the mixture. 
Thus, chain-branching reactions are responsible for the production of pro- 
gressively larger concentrations of chain carriers. Unless the chain- 
branching reactions are balanced by chain-breaking reactions, the reaction 
rate will increase until it becomes sufficiently rapid to be classified as 
an explosion. 

The reaction rate constants are expressed in the standard Arrhenius form 
given by 

= OT^ exp. Y 



KgQ = DCT^^ exp. (FC/T) 

where D, E, F, DC, EC and FC are the constants and are presented in Table 3.1-2. 

The constant, D, for reaction (s) was multiplied by the weighted third body 

efficiency factor, b . 

’ w 

3.2 RESULTS AND DISCUSSION 

3.2. I Effect of Static Temperature 

Chemical kinetic calculations were made as a function of initial static 
temperature for clean air as defined in Table 3.2-1. In all of these calcula- 
tions, the initial static pressure (P = | atm), equivalence ratio (0 = 0.5), 
initial flow velocity (V = 4500 ft/sec) and combustor flow area were held 
constant. 

The results of these studies are presented in Figures 3.2-1 through 3.2-7. 
In Figure 3.2-1, there appears a significant time delay before any appreciable 
temperature rise is observed. This period is referred to as the ignition 
delay period. During this period, little change occurs in hydrogen molecular 
species but large changes in H^O, 0, H, and OH concentrations occur as shown 

in Figures 3.2-4, 3.2-5, and 3.2-6. The two body reactions are relatively 
fast and lead to an overproduction of radicals. At the end of ignition delay 
period, the third body reactions begin to play a dominant role leading to 
recombination and depletion of radicals with the accompanying energy release. 

The time for the process to proceed from the end of ignition delay to an 
equilibrium state is referred to as the reaction time. Due to an Inherent 
difficulty of the computer program in reaching the equilibrium state, the com- 
putations were limited to 3 min of computer time; therefore, the final 
equilibrium state was not reached in most cases. 
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Figure 3.2-1. Effect of Initial Static Temperature on Combustor 
Length 
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Figure 3.2-2. Effect of Initial Static Temperature on Static 
Pressure 
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COMBUSTOR LENGTH 

Figure 3.2.3. Effect of Initial Static Temperature on Chemical 
Efficiency 
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Figure 3.2-5. 


Species Concentration as 
Combustor Length for an 
Temperature 


a Fun< 


Initia 



AIRESEARCH manufacturing COMPANY 

LOS Angeles, Caiiiornia 


UNCLASSIFIED 


s 

70-6319 
Page 3-7 



UNCLASSIFIED 


Hoiijvbj now 



o 

NOIIStfVd S'lOU 


10 

70-6319 
Page 3-8 


AIRESEARCH manufacturing COMPANY 

LosAnc#t«i Cal't(irn.a 


UNCLASSIFIED 


Figure 3.2-6. Species Concentration as a Function of Combustor 
Length for a Static Temperature of I900°R 
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TABLE 3.1-2 (REFERENCE 3-5) 
REACTION-RATE CONSTANTS 


r 

D 

E 

F 

DC 

EC 

FC 

■ 

.23000+14 

.00000 

-.26200+04 

• 

.21592-01 

.28430-00 

..79233+04 

B 

. 10000+15 

.00000 

-.80500+04 

.40779+03 

-.41038-00 

-.86628+04 

E 

.12000+14 

.00000 

-.46300+04 

.18825+01 

.21810-01 

-.91546+03 

B 

. 10000+22 

-.15000+01 

-.51900+05 

.26160+01 

.63998-02 

-.52421+05 

5 

.86000+15* 

.00000 

. 64400+03 

.69905-01 

.28001-00 

.23858+05 

6 

.54000+12 

.00000 

-.12090+05 

.66604-00 

-.46383-01 

-.76370+04 

7 

.29000+15 

.00000 

-.94000+04 

1 

.87164+02 

-.26247-00 

-.88387+04 

8 

.40000+19 

-. 10000+01 

.00000 

.51895-00 

.15772-01 

.51506+05 


. 1 1370+26 

-.25000+01 

-.59400+05 

.78438+03 

-.42620-00 

-.60168+05 

B 

.75000+20 

-. 10000+01 

.00000 

.59569-02 

.27819-00 

.60344+05 

1 

.25000+13 

.0000 

-.19630+05 

.76439+03 

-.38806-00 

-.95775+04 


b = 5.0 X„ + 1.75 X- + 2.5 X„ + 1.00 X.„ + 30.0 X„ . 
w Om AR HaO 


*2.2 + 15 used in calculations reported in Section 6. 

There are a number of quantitative definitions of the ignition delay 
time which can be found in the literature: 

(a) The time required for the concentration of OH radical to reach 
lO”^ mole/liter (Reference 3-2) 

(b) The time at which a point of inflection of atomic hydrogen is 
observed (Reference 3-3) 

(c) The location of the intersection of the slopes of the initial and 
maximum temperature rise (Reference 3-4) 

(d) The time at which a point of inflection of OH radicals is observed 
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(®) The point of 0.1 percent density decrease 

(f) The time it takes for the temperature to increase 5 percent 
of the total temperature rise (Reference 3-5) 

All these assumptions in defining the ignition delay period are reason- 
able and lead to the same conclusion (Reference 3-6). For convenience, 
definition (a) is used in this report. The summarized results are shown in 
Figure 3.2-7. As expected, the ignition delay increases with decreasing 
temperature and does not ignite at all when temperature decreases to I800°R. 
The resultant pressure and chemical efficiency variations are presented in 
Figures 3.2-2 and 3.2-3. The chemical efficiency is an approximation based 
on the percentage of reacted hydrogen molecules by weight. For complete 
combustion, the percentage of unreacted hydrogen molecules increases with 
c^bustion temperature. Therefore, the approximation of chemical efficiency 
will be pessimistic, especially at high combustion temperatures. The approxi- 
mation was used as a calculation convenience and is permissible because only 
qualitative results are being sought. 

TABLE 3.2-1 

MASS FRACTIONS OF REACTING FLOW 
(EQUIVALENCE RATIO =0.5) 


Clean Ai r 

5.05 Percent Vitiation 
Static Temperature 

49.25 Percent Vitiation 
Static Temperature 

2000° R 

2300° R 

2000° R 

2300° R 

H 

1.000-09 

1.000000-09 

1 .0000-09 

1.000000-09 

1.00000-09 

0 

1.0000-09 

1.000000-09 

0.296-07 

1.000000-09 

0.346-07 


0.7569-00 

0.71959-00 

0.71981-00 

0.38336-00 

0.38348-00 

«2 

0. 1440-01 

0. 14634-01 

0. 14634-01 

0. 17194-01 

0. 17194-01 

°2 

0.2287-00 

0.23228-00 

0.23206-00 

0.27292-00 

0.27279-00 

OH 

1.0000-09 

0.345-06 

0.3438-05 

0. 1 16-05 

0.116-04 

«2° 

1 .0000-09 

0.33491-01 

0.33491-01 

0.32653-00 

0.32653-00 

HO2 

1 .0000-09 

1 .000000-09 

1 .0000-09 

1 .000000-09 

0.330-07 

CM 

0 

CM 

X 

1 .0000-09 

1.000000-09 

1 .0000-09 

1 .000000-09 

1 .0000-09 
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3*2«2 Effect of Pressure 


Chemicat kinetics calculations were made for a hydrogen-air mixture at 
a temperature of |900®R for several initial pressures in a constant area duct. 
The ignition delay times obtained from these calculations are presented in 
Figure 3.2-8. At low pressures, the collision rate of molecules Is small due 
to the low density of the mixture, therefore, the chain-branching reactions 
are slow which results in long ignition delays. As the pressure is increased 
more and more molecules collide per second and chain-branching reactions 
become faster and faster resulting in shorter ignition delays. However, when 
the pressure is further increased, the ignition delay time reaches a minimum 
point beyond which it increases with increasing pressure. This phenomenon is 
caused by the following chain-breaking reaction: 

H + O2 + M = HOg + M 

where M represents any kind of molecules, such as H^, 0^, H^O, N^, etc. 
Ordinarily HO^ is a very unstable combination of hydrogen and oxygen because 

more energy is released upon its formation than is required to bring about its 
decomposition. However, if another molecule is present to absorb the excess 
energy released from the formation of HO 2 , the latter may exit as an inter- 
mediate product in the chain reaction. The "third body", M, acts as an energy 
acceptor and stabilizer to this reaction. This reaction serves to break the 
hydrogen -oxygen chain reaction by removing the free hydrogen atom from the 
mixture. As the pressure increases, the concentration of "third body" 
increases and the effect of chain-breaking becomes greater than the effect of 
chain-branching; therefore, it increases the Ignition delay time. 

_Figures^3.2-9, 3.2-10, and 3.2-11 show the temperature, pressure, and 
chemical efficiency distribution along a constant area combustor with various 
inlet static pressures. The concentration of species is given in Figure 3.2-12 
at an inlet pressure of 25 psia. 

The effect of pressure on ignition delay for an initial temperature of-- 
2200 R is presented in Figure 3.2-13. The ignition delay time decreases with 
increasing pressure and no minimum point was observed up to a pressure of 65 
psia. The reaction is slowed down slightly in the range of pressures from 
65 to 100 psia. Figures 3.2-14, 3.2-15, and 3.2-16 show the temperature, 
pressure and chemical efficiency distributions, respectively, along the duct 
with initial inlet temperatures of 2200®R. 

Calculations were also made for other inlet temperatures which are shown 
in Figure 3.2-17. As the static temperature is increased from I900°R to 
2200®R, the pressure at which the minimum delay time occurs moves progressively 
to higher values and the curves become flatter. Thus, at the higher initial 
temperature of 2200®R, ignition delay was not affected as much by static 
pressure as at an initial temperature of |900*R. 
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Figure 3.2-8. Effect of. Initial Static Pressure on Total Ignition 
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Figure 3.2-10. Effect of Initial Static Pressure on Pressure as 
a Function of Combustor Length 




Figure 3.2-11, Effect of Initiel Static Pressure on Chemical Efficiency 
as a Function of Combustor Length 
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significant. Therefore, at low initial static temperatures, vitiation caused 
higher ignition delays. At initial static temperatures above 2200°R the 
difference in ignition delay time for vitiated and clean air were less than 
20 psec and therefore can be considered negligible. 

3.2.5 Effect of Initial Radicals (Carriers) 

In order to see the effect of initial radicals or carriers on the 
chemical reaction as a function of initial static temperature, chemical 
kinetic studies were made for two vitiated mixtures with initially non- 
equilibrium free-radical concentrations (see Table 3.2-2). Results are 
shown in Figures 3.2-23 through 3.2-28 as plots of temperature, chemical 
efficiency, and pressure versus combustor length with different initial 
temperatures for 10 percent and 2 percent of initial carriers. The percent 
of carriers indicates the amount of hydrogen consumed in a constant area 
duct if the initial temperatures were 2200*R. As expected, the presence of 
initial carriers accelerated the chemical reaction by the chain-branching 
reactions of Table 3.1-1, steps 2, 3, and 7. Initial carriers had the effect 
of transposing the region of long ignition delay to lower temperatures. The 
autoignition temperatures (those temperatures below which chemical reaction 
ceased to occur) were lowered to I600®R and I400®R for the vitiated mixtures 
of 2-percent and 10-percent reacted-cha in carriers respectively, as compared 
to I800°R for clean air. Ignition delay could not be determined by the 
previously used criterion of OH concentration because of the presence of 
initial carriers. Results based on the chemical efficiencies of 80 percent 
and 40 percent for both vitiated mixtures and clean air are shown in 
Figure 3.2-29. For a given combustor length. Figure 3.2-29 gives the 
required initial static temperature to produce a given chemical efficiency. 
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Figure 3.2-18. Effect of Equivalence Ratio on Ignition Delay 
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Figure 3.2-22. Effect of Vitiation on Ignition Delay Time 
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Figure 3.2-23. Effect of 10-Percent Carriers on Static Temperature 
as a Function of Combustor Length for Several 
Initiai Static Temperatures 
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Figure 3.2-24. Effect of 10-Percent Carriers on Chemical Efficiency as 
a Function of Combustor Length for Several Initial 
Static Temperatures 
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Figure 3.2-25. Effect of 10-Percent Carriers on Static Pressure as 
a Function of Combustor Length for Several Initial 
Static Temperatures 
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Figure 3.2-27. Effect of 2-Percent Carriers on Chemical Efficiency 
as a Function of Combustor Length for Several 
Initial Static Temperatures 
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Figure 3.2-28. Effect of 2-Percent Carriers on Static Pressure 
9S d Funct t on of Combustor Lgrq th for Severd I 
Initial Static Temperatures 
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Figure 3.2-29. Effect of Percentage of Carriers on Initial Static 
Temperatures as a Function of Required Combustor 
Length and Chemical Efficiency 
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TABLE 3.2-2 

MASS FRACTIONS OF REACTING FLOW WITH INITIAL CARRIERS 
(EquivalenceRatio = 0.5) 



2 Percent Carriers* 

10 Percent Carriers* 

H 

0.63844118-04 

0.30009153-03 

0 

0.34673790-03 

0.16987557-02 

^2 

0.75690000-00 

0.75690000-00 

»2 

0. 14104986-01 

0. 12964130-01 

°2 

0.22605520-00 

0.21569300-00 

OH 

0.88783859-04 

0.44593853-03 

«2° 

0. 18617049-02 

0.9132536-02 

HO 2 

0.58185451-03 

0.28683717-02 

«2°2 

0.12234745-06 

0.64002613-06 


^Percent carriers indicates the amount of hydrogen consumed in a constant-area 
duct if the initial temperatures were 2200**R. 


3.2.6 Effect of the Divergence 

The effect of divergence was calculated for clean air as a function of 
initial static temperature. The results for ah initial temperature of 2500°R 
are presented in Figures 3.2-30, 3.2-31, and 3.2-32 for different values of 
I dA 

^ ^ (change of area ratio per unit length). 


Examination of Figure 3.2-30 indicated that ignition delay increased with 

with increasing divergence^^ and that the reaction proceeded very slowly 

at values of divergence greater than approximately l/lO cm"*. The Increased 
ignition delay with divergence was due to the relative production rates in the 
chain-branching and chain-breaking reactions. At low values of divergence, 
the production rates in the chain-branching reactions were sufficient to offset 
the production rates in the chain-breaking reactions. At large values of 
divergence, the decreasing values of pressure and temperature effectively 
quenched the chain-branching reactions and the production of chain carriers was 
not sufficient to Initiate the highly exothermic recombination reactions. For 
example, to achieve a chemical efficiency of 80 percent within a combustor 
length of 4 cm, the required initial static temperature would be 1670°R, 

I9I0°R and 24I0®R for 10-percent, 2-percent and 0-percent carriers, respectively. 
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Figure 3.2-30, Effect of Divergence on Static Temperatu 
Function of Combustor Length 
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Figure 3.2-31. Effect of Divergence on Chemical Efficiency 
as a Function of Combustor Length 
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Calculations were also made for clean air at an initial static temperature 
of 2200®R. The results presented in Figures 3.2-33, 3.2-34, and 3.2-35 
indicated a trend similar to the higher initial static temperature case wi th 
the reaction essentially quenched for values of divergences greater than about 

6/ 1 00 cm" * . 

Examination of Figure 3.2-33 indicated that ignition would occur in a 

diverging duct when concentration of OH had reached lO”^ moles/liter before 
the temperature had decreased below about |950®R. 

I dA 

In Figure 3.2-36, the ignition delays plotted versus ^ ^ for initial 

static temperatures of 2200®R, 2350°R, and 2500®R indicated that the effect of 
divergence on ignition delay decreased with increasing initial static tempera- 
ture. The effect of increasing divergence was to increase the length required 
for combustion. This increased length was due to (l) the temperature and 
pressure effects on ignition delay, and (2) increased flow velocity. 


j 
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Figure 3.2-34. Effect of Divergence on Chemical Efficiency 
as a Function of Combustor Length 
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Figure 3.2-35. Effect of Divergence on Static Pressure as 
a Function of Combustor Length 
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4. CHEMICAL KINETICS DURING THE IGNITION DELAY PERIOD 

Although more than a dozen chemical reactions are involved in a hydrogen- 
air combustion process, not all reactions are significant at all times. In 
the isothermal induction period the major change is in the buildup of the 
atomic species (H) and (O) and the (OH) radicals. The significant reactions 
are either mildly exothermic or endothermic, and the system is essentially at 
a constant temperature. The production rates in reactions which involve large 
activation energies are insignificant during this isothermal induction period. 
The length of this period is affected both by the initial temperature and by 
the pressure. For the HRE application, the ignition delay must be considerably 
less than one millisecond. Significant production rates occur in the following 
reactions during the induction period: 

H + 0 , *^! - * 0H + 0 k , = 5.64 x lo'° e - (Ref 3-4) (4-1 ) 

= 8.6 X 10® e + *^^® (Ref 3-5) (4-2) 

kj = 1.2 X lo'° e - (Ref 3-4) (4-3) 

H^ + 0H-.!1 ^ -^H^0 + H k^ = 4.2 X loV^ e - ^ (Ref 3-4) (4-4) 

Where the reactance concentration is in units of mole/liter, R is in kilo-cal/ 
mole ®K, T is in °K, and the unit for time is seconds. 

From the above equations the net production rates of (OH), (H) and (O) 
can be readily written: 



= ShW'" ''3l"2ll°l - (*-5) 

,4.7) 
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Observing the numerical solutions for more complete reactions in Section 3, 
the relationship 



remained relatively constant for a given initial condition. Also during the 
ignition delay period, the molecular concentrations did not change materially. 


and the ratio of 


defined as 


order to simplify the solution 
be treated as a constant. 


20 could be considered constant. Therefore, in 
the following quantity denoted by symbol B, will 


k_ fH J 


0 

CM 

ro 



k. 10 J 


TT' 

CM 




Substi tuti 
results in the 

ng these two definitions into Equations (4-5), (4-6), 
following equations. 

and (4-7) 

d[0H] 

dt 

= a, [H] 

- ag 

(4-8) 

d[H] 

dt 

= [oh: 

- a^ [H] 

(4-9) 

d[0] 

dt 

- aj [H] 


(4-10) 

where: a^ s 

N “i " 

+ B) 


"2 ^ 

l°2| '‘4 




hi S (' 

- B) 

(4-11) 

% " 

hits " 

- B) + k^ [M]} 


Equations (4-S) and (4 
oxygen concentration. The 

-9) are coupled and are independent of the 
solution of Equation (4-8) and (4-9) is 

atom! c 

[OH] 

m t 

-S' 

m_t 

"S' 

(4-12) 

[H] ^ 


”l ' SS ”2‘ . b, 

(4-13) 

3/ + 

a, + 
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where b , b and b are constants to be evaluated, and m and are 

12^ ' I o 


-(a^ %) + ^ 


-(a^ ^ %) - + 4a^ (a, - a^) 


The term involving m^ has a large negative value and dissipates rapidly, 

and only mj is of significance. For the range of temperatures and pressures 

important to HRE operation, (900°K < T < I400°K and pressure =■ I ATm) the 
values for aj, a^ and a, are: 


0 = 0.5 


= 1.0 


T = 900®K a I = 2.35 x 10^ (l + B) 


Same as for 0 = 0.5 


a^ = 7.4 X 10' 


1.48 X 10 


a^ = 2.3 X 10^ (l - B) + 3.5 x 10^ Same as for 0 = 0.5 


T = I400®K a I = 4.8 X 10^ (l + B) 


Same as for 0 = 0.5 


a^ = 8.6 X 10 


1.7 X 10 


= 4.8 X 10^ (l - B) + 2.7 x 10^ Same as for 0 = 0.5 

For temperature exceeding I400®K the ignition delay time will be in 
microseconds. 

The value of B for the stated critical conditions is approximately 0.8. 
From the above table, a_ is an order of magnitude larger than a or a . 

^ I 4 

The exponent mj is then approximated 
- - "2 

m, “ 

I =2 + ^ 

- a -a 
I 4 


■ l°2ll''l ^ - K| U - B) . kj [H] 

- |°2||2 S ' - Kj 
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The [oh] concentration is then given by 


[oh] = [OH], e 


|02ll2k,B - 


(4-17) 


Using the definition of ignition delay as the time for the [OH] radical 
to achieve a concentration of I0~^ mole/liter, then 


TTO 7 


° 2 | P‘"®I “ '^2 ’^”^1 


(4-18) 


- p2].|2k,B-k2 LMJ.| 

Equation (4-IS) was plotted in Figure 4.1-1 for pressures of 1/2, I, and 
2 atmospheres for the temperature range between 900^K and I400**K. 

Also shown in Figure 4.1-1 were the shock tube experimental results taken 
from Ref 4-1. The definition of ignition delay used in Ref 4-1 was based on 
the sharp rise in temperature observed by the infra red optical system. If 
this sharp temperature rise occurs near an [OH] concentration of I0~6 moles/ 
liter, then the agreement between theory and experiment is reasonable. 

The oxygen concentration and the molar concentration [M] of the mixture 
can be replaced by pressure in Equation (4-18). 


where: a. 


a^p - a^p 


0.0037Bk. 


'"rw: 


sec gm 


(4-19) 


1.7 X 10"’' k. 


,2 , I0"“ 

T IOT 7 


sec gm 


P in units of gm/cm^ 


T in 
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Equation (4-19) has a ininimuin when the pressure is varied. 


dT 

dp 


^5 - 

■[a^P - a^p2]2 


= 0 


or 


109 Bk. 


P = 


2a. 


= 7100 BTe 


-16.38/RT 


For B = 0.8 the minimum delay time occurs at the following pressures. The 
solutions from numerical integration are also listed for comparison. 


Simplified Numerical 

Solution Solution 

Lz H atm P, atm 

1000 1.45 I 

1140 4.4 2.7 


The disagreement between these two methods may be attributed to the 
ignorance of the third body effects of water vapor in the simplified solution. 

The ignition delay is also affected by the initial [OHj concentration. 

The influence of this factor can be deduced from equation (4-18) as follows: 

T ® (log I0“* - log 

Increasing the initial [OH] concentration wi 1 1 reduce the ignition time. 

A typical initial [OH] concentration Is 10-10 mole/liter at IOOO®K. A tenfold 
increase of this value produces a 25 percent reduction of the delay time while 
a lOO-fold increase results in a 50-percent reduction as calculated below: 


“^10 fold [OH] 
'"[OH], 


= 10 


-6 - (-9) _ 3 

-6 - (-lOj ■ 4 


T 


T 


100 fold [OH] 


[OH]. 


= 10 


-10 


-6 - ^8) _ j_ 

-6 - (- 10 ) “ 2 
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The increasing ignition delay in the second explosion limit is clearly 
illustrated by Figure 4.1-1. For an initial temperature of IOOO°K or less, 
the ignition delay time increases with pressure. This phenomena has been 
known for some time and is due to the formation of the [HO^] radial by 

reaction 2 which is (H + 0^ + M ■* HO 2 + M). This reaction removes the atomic 

hydrogen without the benefit of producing a new active species. Since it is a 
three-body-col I ision process, its rate is proportional to the second power of 
the pressure as compared with the two-body-col I is ion process which is propor- 
tional to the first power of the pressure. The ignition delay time will 
increase with increasing pressure when the rate of reaction 2 increases faster 
than reaction 3. 

The initial rate of [OH] production is primarily given by the molecular 
collision denoted by the reaction 

+ 0^ — 5-^20H 

At time zero, it is assumed that the production rate of [OH] in the above 
reaction balance the [OH] consumption rate of reaction 4. 



The calculated data presented in Figure 4.1-1 were computed based on the 
initial [OH] concentration given by Equation (4-20). 
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5. TWO-DIMENSIONAL EFFECTS IN A DIVERGING COMBUSTOR 

A few exploratory calculations were made In order to Investigate two- 
dimensional effects in supersonic combustion of a premixed gas in a diverging 
cOTbustor. Since this effort was very limited In scope, only a brief qualita- 
tive discussion will be given here. 

Numeriul solutions were obtained using the method of characteristics with 
finite chemical kinetics in an axisymmetric combustor. The computer program 
described in Reference 5-1 and 5-2 was used to obtain these solutions. 

The hot portion of a boundary layer can give beneficial effects in the 
initiation of combustion (assuming that there is fuel in the boundary layer). 
Ignition and combustion in the boundary layer can produce an expansion of the 
boundary stream tube. Compression waves are induced in the main flow where 
the stream tube begins to expand. Farther downstream, when combustion is 
nearly complete in the boundary layer, expansion waves appear in the main flow. 
This situation is depicted schematically in Figure 5.0-1. At the duct center- 
line the pressure waves appear downstream of the start of the compression 
surfaces. For typical AIM Mach numbers this distance is of the order of the 
duct height. 

Compression heating of the mainstream, caused by combustion in the boundary 
tube, was predicted in Reference 5-3 using a quas i -two-dimens iona 1 
model which assumed a uniform pressure across the combustor's height. The 
"overshoot" (See Figure 5.0-1) predicted by the truly two-dimensional solution 
has an effect on combustion which is usually benefical and which is not pre- 
dicted by the quas 1 -two-dimens iona 1 approach. 

However, the beneficial effect of a high initial temperature in the 
boundary layer can be offset. When a diverging portion of the combustor is 
reached, the boundary layer stream tube is the first to be turned (expanded). 
This can result in quenching of combustion in the boundary layer. 

Also, choking can occur in a single stream tube, possibly leading to a 
disturbance which propagates upstream in the combustor and causes an adjust- 
ment in the air flow. Such a disturbance cannot be predicted using the 
existing steady-state model. However, this phenomenon could conceivably occur 
in a real combustor and be a factor in boundary-layer separation. 
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Figure 3.2-15, 


Effect of Initial Static Pressure on 
Pressure as a Function of Combustor 
Length 
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Figure 3.2-17. Effect of Initial Static Pressure and 
Temperature on Ignition Delay 
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3.2.3 Effect of Equivalence Ratio 

Chemical kinetics were calculated as a function of equivalence ratio with 
clean air. Results presented in Figures 3.2-18 through 3.2-21 indicated that 
ignition delays are essentially constant for the range of equivalence ratios 
i nvest i gated. 

3.2.4 Vitiation Effects 


Precombustion of hydrogen with oxygen replenishment as a means of perheated 
air is denoted by the term "vitiated air." The effect of this vitiated air on 
the ignition delay time of the -air reaction was investigated. 

Clean air at several initial total temperatures was precombusted at a 
pressure of 20 atm with varying amounts of hydrogen at room temperature. The 
combustion products were then expanded to a pressure of I atm in order to 
obtain their equilibrium compositions. Note that in the precombustion, oxygen 
was added to preserve the same 0^ mole fraction as in clean air. The percent 

of vitiation or percent of fuel by weight (fuel considered as and 0 ^ with 

the weight of 0_ 12.2 times that of H.) was varied from 0 to 50 percent (see 
Table 3-2-1). ^ 2 

Chemical kinetic calculations were then made starting with the vitiated 
equilibrium species concentrations at I atm pressure. Hydrogen at an equiv- 
alence ratio of 0.5 was allowed to react at constant area with the vitiated 
air. The initial conditions consisted of a static pressure of I atm and a 
velocity of 4500 ft/sec at several initial static temperatures from I850®R to 
2300®R. The results of these calculations are presented in Figure 3.2-22. 
Ignition delays were determined from the time required for the OH concentra- 
tion to reach 10"6 moles/1 i ter. The dotted lines are the estimated values. 

Examination of Figure 3.2-22 indicated that the degree of vitiation had 
a strong effect on ignition delay time at I850°R and 2000®R, but a negligible 
effect at 2300®R. Even small amounts of vitiation at I850®R caused the 
ignition delay time to increase very rapidly and for all practical purposes 
caused chemical reactions to cease. Chemical reactions also ceased at 2000°R 
initial temperature near 25 percent vitiation. 

As the amount of vitiation increased, so did the amount of water vapor 
and OH. The initial concentrations were shown in Table 3.2-1. The water 
vapor acts as a very efficient third body (third bodies act as energy sources 
or sinks in chemical reactions) as shown in the reaction: 

H + Og + M = HO^ + M 

This reaction removed H from the mixture causing a quenching of the 
chain-branching reaction. At low temperature, the chain-breaking reaction 
H + 0^ + M = HO^ + M predominated. At the higher temperatures, the thermal 

dissociation of H^O and chain-branching HO^ + H^ = H^O^ + H then became 
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6. COMPARISON WITH TEST DATA 

Considerable attention has been given to the theoretical chemical kinetics 
of premixed hydrogen-air supersonic combustors. However, little of this work 
has been done for diverging combustors of specified geometry. In this type 
of combustor, conditions can occur in which the expansion waves cause the 
temperature and pressure to decrease from their entrance values with combustor 
length. By contrast, if constant pressure, constant area, or constant Mach 
number is assumed, the temperature never decreases appreciably (unless large 
wall heat losses are considered). 

Because the rates of the elementary chemical reactions in the hydrogen- 
air combustion system are strongly dependent on concentrations and temperature, 
it seems plausible that under some conditions divergence of a supersonic com- 
bustor could lead to significant adverse effects on the chemical kinetics even 
though the initial temperature is above the autoignition temperature. 

Experimental results reported by Cbokson, Flanagan, and Penney (Reference 
2-1 ) provide an example of conditions under which a small change in divergence 
angle had a dramatic effect on combustion. 

Cookson et a1. report that for a combustor inlet air total temperature of 
approximately I900*’K and at Mach 1.98, a flame was initiated in a combustor 
which consisted of a 5-in. -long constant-area section followed by a 25-in. - 
long conical section diverging at a I -deg angle; however, with the same inlet 
conditions no combustion took place within a completely divergent duct, even 
with a divergence angle as small as 0.75 deg. 

Axial injection was used in Cookson's experiments. This eliminates com- 
plicating effects associated with cross-stream injection and makes the experi- 
ments more amenable to one-dimensional analysis. 

In an attempt to explain Cookson's results, chemical kinetic calculations 
were made using two different computer programs, each using a different set of 
chemical kinetic rate constants. The first program used (Zupniks) is described 
in Reference 3-1.; the second (Frey's) in Reference 2-2. The chemical reactions 
and rate constants used in the two programs are listed in Tables 6.0-1 and 

3.1- 2. For initial calculations the higher frequency factor listed in Table 

3.1- 2 was used for reaction 5. 

Because the analytical models used were one-dimensional, there was some 
question as to what initial conditions were appropriate. In the Cookson 
experimental work the fuel and air started out unmixed at the combustor 
entrance. For the chemical kinetic calculations, premixing was assumed. All 
initial parameters were kept constant except the initial temperature, the 
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TABLE 6.0-1 

REACTIONS USED IN THE SHEAR LAYER COMBUSTION PROGRAM 


H 2 O + H 


H 2 O + 0 

H 2 O 2 + H 


H + H + M 
0 0 + M 

OH OH 
OH + H 2 
0 + H + M 
OH -P 0 
OH + H 
OH^OH 


H 2 + HO 2 


HO, + M = H + 0, + M 


H 2 O + M = H + OH + M 


N 2 + M 


N + N M 


NO + H * N + O + h 

Oj + H » 02*^0 + h 

H 2 O 2 + M » OH + OH + h 


N 2 0 


"2 ^ ^2 
NO + 0 


H2O + OH 


N + OH 
NO N 

NO ^ NO 


O 2 + N 

OH + H^O 
H 2 O 2 + H 


H.O + HO 2 = H,0« + OH 


7.5 X 10 
1.9 X 10'* 

1.4 X 10 ’^ 

6.0 X 10 ' ' 

2.0 X 10 '® 
3.2 X 10 " 

1.4 X 10 '® 
1.0666 X 10 '® 

9.6 X 10 '® 
1.59 X lo'®b 

M 

I. 17 X 10*^ 


6 X 10*^ 

I0>« 

8.4 X 10 ’^ 

5.3 X lo‘* 

1.5 X 10 *^ 

10 '^ 

1.8 X 10 ® 

5.6 X 10 ’^ 
3. 18 X 10 ' 



0.0134 

0.0 

0.0 

0.0 

1.0 


B 


0.0 

0.0 

49.2644 

5.0 

0.0 

0. I 
5.19 
0.96671 
24.0 

1.0 
0.0 
0.0 


Reference 


40.0 

3.020 

77.9 

9.0 

1.8 


BacIcMard rate constant • 

AT e RT 


b » 

M 

5.0 + 1.75 + 2. 15 X„ + 

”2 °2 "2 
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combustor geometry, and the equivalence ratio. Cookson reports an air total 
temperature of approximately |870<»K (3366®R) and an air Mach number of 1.98 at 

<:®'"Pe'-ature apparently were about 

1176 K (2N7°R) and 290»K (523»R) respectively. Preliminary calculations were 
made using 0 = 0.50. This is slightly above the highest value U = 0.45) 
given by Cookson. 

6. I RESULTS 

in. 1 ! mixed fuel-air temperature of I I50®K, the calculations 

indicated a chemical efficiency of over 60 percent in the 0.75 degree all- 
di verging combustor. At an initial temperature of I I30®K the chemical effi. 
ciency calculated was less than I percent for the same geometry. 

For the 30-in. duct with a 5-in. constant-area section and a l-deg diver- 
gence angle, the calculated chemical efficiency was 60 percent or greater for 
^peratures of IH50K and above. This indicated that the calculations with 
initial temperatures in the band of approximately III5®-II35®K (2007° -2034° R) 

otcu^Jerl'JrT’r combustion ^ 

having a 5-in. constant-area section, and did not occur 
in the all-diverging combustor. • 

the Choking in the theoretical model at an equivalence ratio of 0.5, 

the in tial temperature had to be more carefully selected. Raising the 
initial temperature just a few degrees above 1 1 |5°K could lead to choking in 
the divergent portion following the constant-area section. 


calculations were made for an equ 
Zupnik's and Frey's programs, with an initia 
two different sets of rate constants listed 
initial conditions are listed in Tabie 6.1-1 
the result found by Cookson— that combustion 
with a constant-area section followed by a I 
combustor having a constant divergence angle 
a constant-area combustor was specified^ the 
This apparently is in agreement with Cookson 
form in the constant-area combustor. 


i valence ratio of 0.33 using both 
1 temperature of l||5°K and the 
in Table 3.1-2 and 6.0-1. The 
. Both programs again predicted 
would occur in the combustor 
-deg divergence, but not in the 
0.75 deg. Furthermore, when 
calculation indicated choking, 
s observation that a shock would 


caseJ^^re^^ZiTL^c-^ temperature prof i les for the two diverging combustor 
(Figure 6 I ?) th<» '^ures 6.1-1 and 6. 1-2. In the case with no combustion 
^Figure 6.1-2) the temperature profiles from the two programs were very 
nearly identical and only one curve is shown. H 3 ms were very 

the computer programs is due to 

in constants which were used rather than the differences 

bonven^eJI^b!*' P'-‘> 9 ram name only for 


the conspicuous features of the case with no combustion is the 

and beatnb" 1?^ concentration which occurs in the results from both programs 
and begins well along the rombustor length, as shown in Figure 6.1-3.^ (A 
c ease in (o) concentration begins farther downstream; see Figure 6.1-4.) 
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Figure 6.1-2. Chemical Kinetics Simulation of Cookson's 
Experiment 


6( 

70-6319 

Page 6- 5 


AIRESEARCH manufacturing COMPANY 

lot Ann«t«s. Caitiorma 


UNCLASSIRED 



UNCLASSIFIED 



Figure 6.1-3. Chemical Kinetic Simulation of Cookson's Experiment 
Mole Fraction OH vs Z 
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TABLE 6.1-1 

INITIAL CONDITIONS ASSUMED FOR COOKSON'S COMBUSTOR 

MASS FRACTIONS 


Species 

Zupnik Program 

Frey Program 

H 

10-’ 

0.0 

0 

10 -’ 

0.0 

■^2 

0.761 

0.761 

«2 

0.00958 

0.00958 

°2 

0.230 

0.230 

OH 

10 -’ 

0.0 

H 2 O 

10"’ 

0.0 

HO 2 

I0-’ 

0.0 

»2°2 

10"’ 

0.0 

NO 

Not Appl i cable 

0.0 

N 

Not Appl i cable 

0.0 

°3 

Not App] icable 

0.0 


Velocity, ft/sec = 4215 

Static pressure, psia = 12.81 

0 s 0.33 


Reaction rates from the Frey program indicate that the reversal in the time- 
derivative of (oh) concentration occurs when the rate of the OH-destroy ing 
reaction 


OH+H2-H2O+H ( 6 .,) 

outweighs the combined rates of the two dominant OH-producing reactions 

Og + H - OH + 0 (6_2) 

H 2 + 0 - OH + H (6_3) 


^ 0 - 
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The (oh) reversal is apparently only a symptom and not the basic cause of the 
quenching of combustion in the analytical study. Exploratory calculations 
showed that when the rate of Reaction (6-1 ) above was suppressed by lowering 
its rate constant, the destruction of (OH) was accomplished by the reaction 

OH + OH - H^O + 0 (6-4) 

Similarly, when in a further calculation Reactions (6-1 ) and (6-4) above were 
both suppressed, two more reactions took their place to cause net destruction 
of (OH). 

It is believed that the quenching observed in the computer programs is 
due to a decrease in the rates of all the reactions, which is caused simply by 
the lowering of static pressure and temperature associated with expansion of 
the gas. The OH reversal seems to be a part of the adjustment of the chemical 
process to the slower rates. Changing the rate constants of the OH-destroy ing 
reactions changes the overall process rate, but apparently does not change the 
mechanism. 

The pressure profile reported by Cookson for 0 = 0.33 and those computed 
by the two chemical kinetic programs for the geometry in which combustion 
occurred are shown in Figure 6.1-5. Computed profile forms in disagreement 
with the experimental profile is not surprising; the heat release was spread 
out due to finite-rate mixing in the experiments, but in the theoretical pre- 
mixed model the heat release occurred over a short axial distance. Agreement 
in the constant-area section would probably be improved by considering the 
effects of mixing and of wall friction. 

As previously mentioned, a high frequency factor was assigned to the rate 
constant for reaction 5 in using the Zupnik program. This frequency factor 
was about 2.5 times greater than that in the original reference (see Table 
3.1-2). Decreasing the frequency factor of this reaction has the effect of 
increasing the overall heat release rate. When the lower frequency factor in 
Table 3.1-2 was used, the heat release was accelerated enough to cause choking 
In the combustor having a constant-area section although the calculation for 
the conical combustor still indicated no heat release. Thus for a marginal 
case the results can be very sensitive to the rate of the H + 0^ + M reaction. 

This rate constant Is believed to be known within a factor of two (Reference 

6 - 1 ). 


The temperature chosen for obtaining qualitative agreement with Cookson 's 
results was about 1I0°R, or 5 percent lower than the unmixed air inlet static 
temperature. This is not considered a serious discrepancy in view of the 
assumption of premixing in the model and the uncertainties in the chemical 
kinetic rate constants. 

It is believed that the qualitative agreement between theory and experi- 
ment In this case Is sufficient to indicate that the 0.75-deg divergence of 
the Mmbustor used by Cookson et al. may have had a serious effect on the 
chemical kinetics, leading to the observation of little or no burning Inside 
the combustor. 
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Figure 6, 1-5. Pressure Profiles for Cookson Combustor 
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6.2 


CHEMICAL KINETIC SIMULATION OF TWO-DIMENSIONAL COMBUSTOR TESTS 


Combustion tests were conducted with a two-dimensional, constant-width 

^ length of 29 In., .0 inl.t cross section of 0.6 by 6 in. and 
an exit-to-iniet area ratio of 2.6. The tests were conducted at various test 
temperatures and injector geometries. Fuel was injected in two stages; the 

'"hi second in a diverging area as shown in Figure 

6.2-1. High chemical efficiencies were obtained in the first stage (about 90 
percent) whi le second-stage efficiencies were lower (about 60 to 70 percent). 

A study was made to determine if the second-stage combustion was chemical- 
kinetics limyed or mixing limited. The Zupnik computer program (Reference 

anVthf Fr^! ^ second-s tap combustor p investigate the chemical kinetics, 
and the Frey computer program (Reference 2-2) was used to investigate both 
chemical kinetics and mixing in the first- and second-stage combustors. 

studied was at a test total temperature of 3000° R with 
0 = 0.68 injected in the second stage. Entrance conditions to the second staqe 

chemisJir ® Of’e-dimensional data reduction program using equilibrium 

chemistry. In this case, second-stage entrance conditions were determined 
assuming complete combustion in the first stage. 

The Zupnik computer program (Reference 3-1 ) was used to evaluate second- 
spge combustion. The calculated air static temperature (from the data reduc- 
tion program) was averaged with the I500°R H 2 temperature and resulted in a 

2I00°R mixed fuel-air temperature at the second stage entrance. Results in 
igures 6.2-2 through 6.2-4 indicated no second stage combustion. The OH 

pneentption reached a maximum of 1.47 x I0“® moles/liter or about |00 times 
lower than that used to indicate ignition delay. Figure 6.2-4. 

to thlT t^TonooJ temperature were run in a manner similar 

equivalence ratio of 0.25 was injected in the 
stao^ * ^ respectively, were injected in the second 

It ^Va ®r® 6.2-5 and 6.2-6. In both cases, the fuel 

started burning about I -in. from the second-stage entrance. 

The static temperature distribution computed for complete mixing was 
higher than that calculated from the data reduction program. This higher 
temperatpe arises because more fuel reacts as shown by the chemical efficien- 
cies in Figure 6.2-7. Also, cooling losses are not accounted for. The maxi- 
mum reaction rate peurred (Figure 6.2-7) at a point about 2-in. downstream 
of s®cpd-spge injection, a chemical efficiency of 80 percent was attained 

injection. Several species are plotted in Figures 6.2-8 and 6.2-9 
for the stoipimetric and fuel lean cases, respectively. In both cases there 
was a rapid increase in the concentrations of H, 0, OH and HO^. About 1.5 to 

2.5 in. downstream of injection the production rate of OH, HO^ and H began to 
decrease. 


Therefore with the assumption of complete and instantaneous mixing and 
second-stage air-fuel injection disturbance, the second stage fuel reacted 


no 
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Figure 6.2-1, Combustor Configuration Schemati 
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rapidly at 3800°R total temperature and did not react at all at 3000°R total 
temperature. At 3800®R about 80 percent of all fuel injected reacted which 
was higher than that calculated from the one-dimensional data reduction program. 
A more realistic study at both temperatures was made with another computer 
program. Measured static pressures, cooling, mixing, and friction losses were 
included in the model. The results of that study follow. 

Two of the two-dimensional combustor tests were simulated in the Shear 
Layer Combustion Computer Program (Reference 2-2) developed by Dynamic Science. 
This program permits the analysis of complex chemical kinetic systems with 
arbitrary gas phase reactions and species using a one-dimensional (area or 
pressure defi ned) stream-tube model. Unique features of this program are (l) 
its ability to integrate the differential equations governing the chemical 
kinetic system with complete numerical stability using an implicit integration 
technique (Reference 6-2), and (2) Its ability to add (or extract) mass momen- 
tum, energy, and species at any position along the stream tube. 

The simulation of the tests was based on the assumption of two limiting 
conditions. The first of these conditions assumed that the combustor effi- 
ciencies ulculated from test data represented chemical efficiencies (lOO 
percent mixing efficiency). The second assumed that the combustor efficiencies 
represented the mixing efficiencies. 

The several mixing models used with the two limiting conditions were 
implemented by the concept of an inert hydrogen molecule denoted by the symbol 
HZ^. This inert hydrogen molecule participated in the flow dynamics, but 

was not permitted to dissociate or react with other species. Mixing was done 
by replacing inert hydrogen (HZ^) with active hydrogen (H^) using the species- 

addition feature of the computer program. 

The three basic mixing models are shown in Figure 6.2-10. Model I pro- 
duced early mixing. Model 2 produced a maximum mixing rate at the mid-point 
of the mixing length, and Model 3 produced late mixing. Further variation 
within each mixing model was obtained by varying both the mixing length and 
the amount of hydrogen mixed at the injectors. 

The two tests chosen for this simulation are listed In Table 6.2-1 along 
with their nominal test parameters. 

TABLE 6.2-1 


Time (sec) 

9 

Tj (°R) 

T (“R) 

Combustor Efficiency 


Run 

238 242 

520.30 682.55 

0.6049 1.0267 

2838 3575 \ 

J Station 7.3 
1897 2158 ) 

60 percent 60 percent 
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The chemical kinetics calculations were started at Station 7.3^ just ahead of 
the first stage injectors, using equilibrium concentrations of the species. 

The energy and momentum defects due to heat transfer and wall friction were 
distributed linearly along the combustor, while the energy, momentum, and mass 
addition due to fuel injection were input at the appropriate injectors. The 
pressure-defined stream tube model was employed using measured static pres- 
sures from the two tests (see Figure 6.2-1 i). 

The results obtained from imposing the two limiting conditions on the 
simulation of run 238 in the Shear Layer Combustion Computer Program are pre- 
sented in Figures 6.2-12 through 6.2-14, where the areas, mixing efficiencies, 
and chemical efficiencies were plotted as a function of axial position. In 
all cases 100 percent mixing efficiency (T]^) was assumed in the first stage 

combustor using mixing model 2 with 50 percent mixing at the injectors and 
complete mixing at station 12. Examination of the calculated areas in 
Figure 6.2-12 indicated a reasonable blockage (30 percent) in the constant- 
area portion and reattachment of the flow near the entrance of the divergent 
portion of the first stage combustor. The first stage results demonstrated 
chemical efficiencies of 95 percent at station 13 and nearly 100 percent 
beyond station 17 (see Figure 6.2-14). 

The first limiting condition (lOO percent overall mixing efficiency) was 
imposed on the second-stage combustor using both a 50-percent local mixing 
efficiency at the injectors with mixing Model 2 and a 0 percent local mixing 
efficiency at the injectors with mixing Model 3. An overall mixing efficiency 
of 100 percent was set at station 29 with both mixing models — the former 
simulating early mixing; the latter, late mixing. 
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□ RUN 242 - TIME 682.55 
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Figure 6.2-12. Geometric and Calculated Areas vs Axial Position 
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RUN 238 0= 0.6049 TIME = 520.30 
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Figure 6.2-13. Mixing Efficiency vs Axial Position - 
Run 238 
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RUN 238 <f> - 0.6049 TIME = 520.30 

100^ MIXING EFFICIENCY 
MIXING LENGTHS 

+ 1ST STAGE 3.5 IN. 1 =1.00 

c 

+ 2ND STAGE 7.82 IN. = 0.92 

02ND STAGE 8.2 IN. 1 = 0.82 

c 

60^ MIXING EFFICIENCY 

2ND STAGE MIXING LENGTHS 

A 4.0 IN. *7 = 0.97 
c 

□ 8.2 IN. 1? = 0.96 

c 
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Figure 6.2-14. Chemical Efficiency vs Axial Position 
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Examinations of the calculated areas (symbols + and©) in Figure 6*2-12 
indicated a reasonable blockage (l5 percent) at the second*>stage injectors, 
however both mixing models eventual ly required flow areas greater than the 
geometric areas. Since such flows would be impossible to attain with the 
measured pressures from run 238, the assumption that the combustor efficiency 
equalled the chemical efficiency with 100-percent mixing was in error. This 
argument was further substantiatecl from the chemical efficiencies in Figure 
6.2-14 where values were 92 percent for early mixing and 82 percent for late 
mixing, compared with the value of 60 percent required to match the test data. 

Therefore, the other limiting condition (combustor efficiency = mixing 
efficiency) was imposed on the second-stage combustor. A 50-percent local 
mixing efficiency at the injectors with mixing Model I (early mixing) was 
assumed. An overall mixing efficiency of 60 percent was set at 4.0- and 8.2- 
in. downstream of the Injectors (second stage) as shown in Figure 6.2-13. 

Again, as shown in Figure 6.2-12, the calculated flow areas indicated a 
reasonable blockage in the vicinity of the second-stage injectors, with the 
blockage gradually decreasing to zero as the flow approached the combustor 
exit. The mixing length of 4 in. produced flow areas which closely matched 
the geometric areas and resulted in a chemical efficiency of 97 percent. 

The two limiting conditions were imposed on the simulation of run 242 
using the same mixing models as in run 238 with the exception that time did 
not permit the late mixing case to be run. The flow areas, mixing efficien- 
cies, and chemical efficiencies were plotted as a function of axial distance 
and presented Ln Figures 6.2-15 through 6.2-17. Examination of the areas in 
Figure 6.2-15 again indicated that the first limiting condition (chemical 
efficiency = combustor efficiency, 100 percent mixing) was in error and that 
the flow areas resulting from the second limiting condition more closely 
matched the test; conditions. The fluctuating flow areas near the exit of the 
combustor (Figure 6.2-15) were due to pressure fluctuations (see Figure 6.2-10). 

While it appeared that the second limiting condition, i.e., incomplete 
mixing produced the correct flowfields, the true condition could be found 
somewhere between the two limiting conditions, quite close to the second. If 
the chemical efficiency is denoted by T1 combustor efficiency by T1 and, mixing 

efficiency by Tl^, the relationship. 


T1 


10™ • “n. 


’m 


(6-5) 


could be used to iterate on a mixing efficiency to produce flows which would 
more closely match the test results. For example, the next value of mixing 
efficiency in the iteration for run 238 would be obtained from 




(n + I ) 



0.60 

0.97 


0.62 
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Figure 6.2-17. Chemical Efficiency vs Axial Position 
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The results from the above test simulations support the conclusion that 
the low combustor efficiencies experienced in the two-dimensional combustor 
tests were due to poor mixing in the second stage. 

The detailed results from the chemical -kinetic simulation of runs 238 
and 242 are presented in Appendixes I and 2. 
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7. AIM ENGINE 


7. I ANALYSIS 

In the AIM combustor^ fuel is injected from several locations as 


shown below: . . 

.lb 4 2a 



The first stage injectors (la, and lb) are located in a slightly diverging 
section; the second stage injectors (ic, 4 and 2a, 2c) at greater divergence. 
Fuel flow requirements in the first stage are set to prevent inlet unstart 
and thermal choking. The second stage fuel flow schedule was to be determined 
from cycle analysis and mixing studies. 

A simulation of both stages was made in a two-dimensional combustor 

_ Testing demonstrated high chemical efficiency (near 90 percent) for 
the first stage injectors. The second stage injectors, which were located 
in a diverging area, gave much lower efficiencies (hear 60 to 70 percent). 

A study was made of the present AIM combustor to see if these low efficiencies 
were inherent in the design due to chemical kinetics. 

It was assumed that combustion would be temperature-controlled. Thus 
complete mixing was assumed. The Zupnik computer program was used to de- 
termine whether the second stage combustor could have a high chemical 
efficiency. 

The first case studied was Mach 7 total temperature at 1000 psia total 
pressure to simulate test conditions for future AIM tests. It was assumed 
that an equivalence ratio of 0.4 would be injected from the first stage 
and that the second stage would begin midway between injectors 2a and 2c. 

A chemical efficiency for the first stage was based on two-dimensional com- 
bustor test results which indicated that at least 80 percent chemical ef- 
ficiency should be obtained in the first stage combustor at an equivalence 
ratio of 0.4. 

The engine cycle program was used to determine second-stage entrance con- 
ditions. This program uses equilibrium chemistry. Using the output static 






temperature from the engine cycle program and assuming 900°R fuel at an 

equivalence ratio of 0.6 gave a mixed fuel-air temperature of 3000°R. This 
temperature included second stage cooling losses. Using this mixed tem- 


perature at 19 psia and 4800 ft/sec as the first stage exit conditions, the 
static pressure and temperature profiles for the second stage combustor were 
calculated in the Zupnik chemical kinetics program as shown in Figure 7.1-1 
and 7.1-2. Examination of these profiles indicated a very short ignition 
delay period at mixed temperature of 3000®R. At a distance of 2.7 cm down- 
stream of the second-stage entrance, the flow thermally choked after 85 percent 
of all fuel injected in both stages (0= I.O) had reacted. Instantaneous 
complete mixing was assumed at the second stage entrance. Actually, mixing 
will take some finite length with mixing and reaction occurring simultaneously 
over a portion of the combustor. Therefore, whether choking would actually 
occur cannot be determined using this model. However, with the assumption of 
complete and instantaneous mixing it does appear that at Mach 7 conditions 
there will be no chemical kinetics problem in the second stage combustor of 
the AIM engine. When the first stage chemical efficiency was reduced to 
50 percent, which produced a mixed temperature of 2400®R, the fuel reacted 
less rapidly as indicated in Figure 7.1-2. In this case, a chemical ef- 
ficiency of 83 percent was attained at a distance of 8.3 cm from the second 
stage entrance. Since the Mach number was about 1.6 at this point, it is 

improbable that thermal choking would have occurred had the calculation been 
continued. 


Combustion occurred when species H and OH produced from first stage com- 
busion were removed at the second-stage entrance (see Figure 7.1-2). There 
is, however, a large ignition delay with about 33 percent of the second-stage 
fuel reacted 19.2 cm from the entrance. This case is not very realistic, 
but was conducted to demonstrate that, at Mach 7, there is no combustion' 
problem due to temperature alone. 

An idea of how the mole fractions of several species varied with combustor 
iength is presented in Figures 7.1-3 and 7.1-4. At a mixed temperature of 
3000 R (see Figure 7.1-3) the concentrations of H, 0, OH and HO 2 increased 

rapidly, but after I cm the concentrations of H, 0, and HO 2 decreased. The 
buildup of OH and HgO became less rapid after I cm, and OH started to de- 
decrease Just before thermal choking occurred. Similar trends were observed 
at a mixed temperature of 2400®R (see Figure 7.1-4). It appeared that the 
production rate of HgO started to decrease at 6 cm of combustor length and 

that the reaction of Hg proceeded very slowly beyond 7 cm. In the AIM com- 
bustor there is about 5-times as much length compared with the two-dimensional 
combustor so slower reaction should not be a problem. The variation of second 
stage chemical efficiencies are shown in Figure 7.1-5. The chemical efficien- 
cies were approximated from the change in weight fraction of H^. This is not 

exact since even if the chemical reaction were to go to equi 1 i br ium, some 

would remain. It does indicate, though, a rapid increase tn chemical ef- 
ficiency near 80 percent with a slow rise thereafter. 
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The second case simulated Mach 6 flight conditions. The manner in which 
the study was conducted was the same as at Mach 7. This time an equivalence 
ratio of 0.2 was assumed for the first stage with a chemical efficiency of 80 
percent at a point midway between injectors 2a and 2c. The second-stage 
equivalence ratio was 0.8, and several fuel temperatures from 900°R to I500°R 
were used. The results obtained from the finite kinetics program are shown in 
Figures 7.1-6 and 7.1-7. At a fuel temperature of 900®R, the mixed fuel -air 
temperature was 2l50®R and no temperature rise was observed; at I500°R fuel 
temperature or 2270®R mixed fuel-air temperature the fuel burned rapidly. 

A mixed temperature between these two was run (2200®R) and in this case the 
fuel also reacted. Thus using species concentrations based on equilibrium 
chemistry for second-stage entrance conditions showed that a change of 50°R 
mixture temperature would produce a go or no-go situation. When the H and OH 
produced from the first stage combustion were removed at a mixed temperature 
of 2200 R, no temperature rise was observed, as shown in Figure 7.1-7. This 
was not the case at Mach 7 at a mixed temperature of 2400°R. At a mixed tem- 
perature of 2270 R, about 70 percent of all fuel injected (0 = I.O) reacted 
before thermal choking occurred, while at 2200®R mixed temperature 60 percent 
reacted before computer cut-off time was reached (flow did not choke). The 
chemical efficiencies for these cases are shown in Figure 7.1-8. When the 
fuel reacted it did so rapidly. The chemical efficiencies continued to rise 
at both mixed temperatures* 

Several species were plotted in Figures 7.1-9 to 7.1-12 for cases with 
fuel reacting and without fuel reacting. The OH concentration for a mixed 

temperature of 2I50°R reached a maximum of 2 x lo”^ moles/liter at about 6 cm 
dovmstream of the second stage injection as shown in Figure 7.I-II. At this 
point the static temperature had dropped to about I950®R with a static pressure 
of 20 psia. The parametric studies seemed to indicate that once the static 
temperature became less than about 2000°R in a divergent duct and the OH con- 
centration had not reached 10 moles/liter by then, combustion did not occur. 
Thus, results from chemical kinetics simulation of the AIM configuration were 
consistent with the parametric study in Section, 3 of this report. For the 

cases with combustion(Figures 7.1-9 and 7.1-10) an OH concentration of 10 ^ 
moles/liter was always reached before the static temperature was reduced to 
2000 R. As shown in Figure 7. 1-12, the 50®R difference in static temperature 
between the go and no-go cases caused an increased production rate of OH start- 
ing at about 0.05 cm from the second stage entrance. Thus the OH concentration 

was able to reach 10 moles/liter before the temperature became too low due to 
d i vergence. 


For a mixed temperature of 2200°R, an OH concentration of lO'*^ moles/liter 
was reached in about 3 cm. This corresponded to an ignition delav time of 
about 25 p.sec. By this time the temperature had dropped to 2100®R. if one 
had assumed expansion with friction at a mixed temperature of 2I50°R, the 
variation of static temperature with time would have been as shown in Figure 
7.1-13. With a friction coefficient of 0.006, at 25 p,sec, the static temper- 
ature would be close to 2l00°R. Thus, had friction been included in the model, 
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Figure 7.1-8, Second Stage Chemical Efficiency at Mach 6 - AIM Combustor (U) 
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Figure 7.1-10, Species Concentration at Mach 6 - AIM Combustor (U) 
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Mixed Temperatures at Mach 6 - AIM Combustor (U) 
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combustion at a mixed temperature of 2I50^R would probably have occurred. 

Also, the model did not include any temperature rise from fuel disturbance 
which would have helped initiate combustion. 

The above study indicated that second stage combustion was very sensi- 
tive to the first-stage chemical efficiency, fuel temperature, and skin friction. 

Second-stage entrance conditions (including the concentrations of all 
major species) were determined from an engine cycle program which used equilib- 
rium chemistry. Therefore, a study was made to determine whether application 
of a chemical kinetics solution in both stages would affect the combustion 
characteristics of the second stage. The results showed that i( would. 

The first-stage entrance static temperature at Mach 6 is about I300°R. 

At this temperature the finite kinetics program would indicate no combustion. 

To simulate the air shock due to fuel injection, an air static temperature of 
2000®R was chosen as the first stage entrance condition. An equivalence ratio 
of 0.2 was assumed and combustion at a constant area was specified to insure 
that the fuel would react. In this manner it was hoped that at the same amount 
of fuel reacted (in a constant area), static temperature, pressure and velocity 
would match that from the engine cycle computer program. Thus the only dif- 
ference in second stage entrance conditions (if any) would be the concentration 
of the individual species. 

The first stage was defined to be halfway between Injectors Ic and 4 at 
a chemical efficiency of 50 percent. With 900®R H 2 fuel, the mixed fuel air 

temperature was I800°R. As shown in Figure 7.1-14, reaction took place al- 
most immediately, and the flow choked about I cm downstream of injection. 

In this case the initial concentration of OH at the entrance of the second 

stage combustor was 3 x lO”^ moles/liter, or 30-times that concentration used 
to define the ignition delay time. This OH concentration was also 100 times 
that obtained from the cycle program output which used equilibrium chemistry 
(Figures 7.1-9 to 7.1-12). As shown in Table 7.1-1, the values of carriers 
0, H, HO 2 and H 2 O 2 also changed significantly. 


Therefore with finite rate kinetic chemistry, second stage combustion 
occurred at a mixed temperature of I800°R, whereas with equilibrium chemistry 
a mixed temperature of 2200®R was required for combustion. This indicated 
that at Mach 6, second stage combustion would not be chemical kinetics -1 imi ted. 
Thus at both Mach 6 and 7 there should be no combustion problem if complete 
mixing is assured. 

Since combustion was shown to be very sensitive to initial specie con- 
centration, an investigation was made varying the percent of fuel reacted 
in the first stage in order to determine the effect on combustor length. To 
simplify results the first stage was run at constant area and the second stage 
at a constant divergence. An equivalence ratio of 0.2 was used in the first 
stage at 2000°R, 20 psia and 4400 ft/sec mixed conditions. An equivalence 
ratio of 0.8 was used for the second stage. The effect of varying the per- 
centage of first stage fuel reacted on combustor length is shown in Figure 7.1-15. 
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TABLE 7.1-1 

COMPARISON OF INITIAL MASS FRACTIONS USED IN AIM 
STUDY AT SECOND STAGE ENTRANCE* (u) 


Specie 

Equi librium 
Chemistry 1st Stage 
T)C| = 0.5 

Finite Kinetics 
Chemistry 1st Stage 
T)C| = 0.5 

Equilibrium Chemistry 
1st Stage Tjc. = 0.8 
Used in Figs. 7.1-6 
Thru 7.1-12 

H 

0 

1 

0 

V 

0.24267 X lO’^ 

1 X 10 "^ 

0 

A 

0 

1 

0 

0.35479 X 10 "^ 

1 X 10 “^ 


0 . 74674 

0.746219 

0.7471 1 

«2 

0.02550 

0.025491 

0.02383 

°2 

0.20243 

0. 190273 

0. 18854 

OH 

4.02 X lO"^ 

0.95006 X 10 "^ 

0. 1 1 X 10 ’^ 

H 2 O 

0.02533 

0.0188078 

0.04051 

HO 2 

0 

1 

0 

V 

0.0144648 

1 X lO"’ 

^2°2 

< io-'° 

0.2592 X 10 "^ 

1 X lO’’ 


* 0 ,= 0 . 2 , = 0.8 
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In a length of 10 cm, 2 percent of the first stage fuel reacts while it takes 
about 16 cm for 97.8 percent of the fuel to react. Mixing the temperature of 
output air with 540®R fuel shows that about 25 percent of the first stage fuel 
must react in order for combustion to occur in the second stage at a divergence 

(dA/A|dx) of 0.03 cm * . No combustion occurred in the second stage for values 

of less than 10 percent of fuel reacted in the first stage. However, the 
overall combustor length does not necessarily decrease as the percentage of first- 
stage fuel reacted increases. The reason is that first stage combustor length 
increases much more rapidly than the second stage length. This is more clearly 
shown in Figure 7.1-16. In fact a minimum occurs near 60 percent of first stage 
fuel reacted. The overall length is shown at a constant Hg mole fraction of 

0.15 or at a constant overall chemical efficiency of about 53 percent. This 
efficiency was the highest reached without choking the second stage with 97.8 
percent of fuel reacted in the first stage. A similar minimum occurred when 

the divergence was doubled to 0.06 cm"* (Figure 7.1-17). This indicated that 
making the first stage of the combustor less efficient might decrease the 
overall combustor length. 

In Summary, 

(1) Equilibrium chemistry and finite kinetics chemistry in the 
first-stage combustor did not predict the same combustion 
characteristics for the second stage of the AIM combustor. 

(2) Finite-rate chemistry assuming complete mixing indicated 
no second stage combustion problems at both Mach 6 and 7. 

(3) Shorter combustor lengths might occur with less efficient 
first stages. 

7.2 AIM KINETICS AT MACH 6 

Model T inlet test data indicate that at Mach 6 the static temperature 
for undisturbed flow at the first injecti.on station of the AIM will be 
approximately |290°R, the local Mach number being about 2.8. This static 
temperature is too low for autoignition. Some exploratory chemical kinetic 
calculations were made, using the Shear Layer One-dimensional Combustion 
Program, (Reference 2-2 ) to ascertain what strength of flow disturbance might 
be needed at the first injection station in order to initiate combustion. A 
further calculation was made for the second stage with no flow disturbance 
other than that from simple mass and energy addition (no boundary layer sep- 
aration). 

For ail these calculations certain constant conditions were assumed. 

The first-stage total equivalence ratio was set at 0.2, the fuel being 
equally split between the two first stage injectors (la and lb). The second- 
stage total equivalence ratio was 0.8, split equally between Tnjectors 2c 
and 2a. The igniters, about I in. downstream of the first-stage injectors, 
were assumed to contribute 0.0413 Ib/sec of hydrogen-oxygen combustion products 
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Figure 7.1-16. Combustor Length vs Amount of First Stage Fuel 
Reaction (u) 
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Figure 7.1-17. Combustor Length vs Amount of First Stage Fuel 
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at an overall igniter equivalence ratio of L6. The igniters* combustion 
efficiency was taken to be 23 percent, giving a mixed igniter exit gas tem- 
perature of 3050®R. The composition of the igniter gases was taken from quasi- 
equilibrium calculations using the above chemical efficiency. 

A uniform frictional force of 18.7 Ibf/in. along the combustor’s length 
was assumed. Estimated local heat losses were used; these varied between 109 
Btu/in. sec. and 234 Btu/in. sec. and were specified as a function of axial 
distance. 

The fuel temperature was assumed to be I440®R, and the total air flow was 
taken to be 34.6 Ib/sec. The fuel and the igniter gases were assumed to be 
injected normal to the engine axis, and to mix completely and instantaneously 
with the main flow. The fuel was injected over a length of 0.2 inches at each 
i njector. 

The disturbances assumed for the flow ahead of the first Injection station 
corresponded roughly to those produced by two-dimensional flow over a wedge. 
Thus the strength of the disturbance can be described by specifying a wedge 
angle. 

Resul ts 


Assuming a 15-deg wedge angle ahead of the first Injection station gave 
a Mach number of 2.1, a static pressure of 47 psia, and a static temperature 
of I740®R just upstream of the first injector. An effective flow area reduction 
of 47 percent was used in order to satisfy continuity. In the first calculation 
using these conditions, the effective flow area was kept constant for a distance 
of 2 in., then linearly increased to the actual duct area at. a point 8 in. down- 
stream of the first injector. Fuel Injection and igniter gas injection pro- 
duced static temperature rises in the constant-area section (see Figure 7.2-1), 
but the maximum temperature was less than I800®R. Ignition did not occur. 

The computed pressure profile is shown in Figure 7.2-2, and Figure 7.2-3 shows 
the concentration of OH. 

A second calculation was made using the same initial conditions but with 
the assumption that the pressure along the combustor was constant. In this 
case fuel injection produced a slight drop in static temperature, and the 
effect of the igniter on static temperature was negligible (see Figure 7.2-4). 
Again, ignition did not occur. 

Increasing the wedge angle to 20 deg gave an initial static temperature 
of I940°R, a pressure of 61.7 psia, and an effective area reduction of 51 
percent. A chemical kinetic calculation using constant pressure indicated 
complete combustion in the first stage. The temperature profile is shown in 
Figure 7.2-5. The area expansion of the gas stream at the equivalence ratio 
of 0.2 (see Figure 7.2-6) was not sufficient to fill the actual cross-sectional 
area of the AIM duct. Consequently, the calculation was re-started at a point 
near completion of first stage combustion (Z » 45; see Figure, 7.2-5). From 
this point on, an assumed effective flow area (rather than constant pressure) 
was specified. The initial area was that determined by the first stage constant 
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Figure 7.2-2. Computed Pressure Profile(U) 
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Figure 7.2-3. Concentration of 
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Figure 7.2-4. Effect of Igniter on Static Temperature(u) 
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Figure 7.2-5. Temperature Profile(u) 
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pressure calculation; the effective area was then increased linearly until it 
matched the actual duct area just downstream of the first injector of the 
second stage. At the entrance of the second stage, the computed gas temperatu 
was over 2300°R. The calculation indicated that the fuel from the first 
injector of the second stage would burn in a very short distance and cause 
choking of the flow. 

7.3 DISCUSSION 

Time limitations prevented the investigation of any mixing model other 
than instantaneous mixing for the AIM. With respect to the mixing model, the 
calculation must be considered optimistic. Also, it is not known whether any- 
thing approximating a constant-pressure profile can be sustained over a length 
of 4.5 in. in the AIM's diverging first stage geometry. 

The strength of the flow disturbance (pressure ratio of almost 4 to l) 
required to produce combustion in the theoretical calculations is probably 
excessive; it is unlikely that this pressure ratio can be induced in the AIM 
(at some equivalence ratio) without unstarting the inlet. However, it may be 
possible to initiate combustion with smaller disturbances because of localized 
igniter effects. 

Treating fuel and igniter gas injections one-dimensional I y undoubtedly 
neglects some of the most beneficial effects of these processes. The igniter 
gases, for instance, are known to be poorly mixed; the igniters* design is in 
fact based on film cooling by unmixed oxygen. The core of the igniter streams 
therefore contains zones which have higher temperatures, and correspondingly 
higher concentrations of labile species, than indicated by the fully-mixed 
igni ter model . 

The effectiveness of the igniters is also considerably increased because 
the hot gases mix only with part of the mainstream. The assumption used in 
the calculations, namely that mainstream gas is mixed with igniter gases in 
a mass ratio greater than 800 to I, leads to the unrealistic result that the 
igniters* effect on local static temperature is negligible. For Mach 5 test 
conditions the igniter used in the two-dimensional combustor was effective 
with a static temperature (before boundary layer separation) of approximately 
I000°R. The undisturbed static temperature in the AIM is I290®R; hence, it 
seems possible that combustion can be initiated in the AIM by use of the 
igniters even if the flow disturbance is small. Also, the prediction of choking 
at the second injection station might be negated by use of a reasonable finite- 
mixing model . 

The conclusion from this AIM study is that with instantaneous mixing and 
for the conditions assumed in the 20 deg wedge angle case, the second stage of 
the AIM would not be kinetically limited at the Mach 6 test condition. 
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8. DISCUSSION AND RECOMMENDATIONS 


8.1 DISCUSSION 


The analytical studies have indicated that chemical kinetics problems 
would not be present in the AIM engine at Mach 6 and 7 flight conditions. 
Inadequate mixing was considered to be the major contributing factor for the 
low combustor efficiencies in the two-dimensional combustor tests. In order 
to substantiate these theoretical findings, it is recommended that future 
research programs be directed to better understanding the mixing process, par- 
ticularly in a confined, favorable pressure gradient system. 


As discussed in the fifth oral review. Reference 8-1, the mixing process 
in a supersonic combustor may be classified in three distinct phases; (l) jet 
penetration, (2) turbulent mixing, and (3) molecular diffusion. Among these, 
the turbulent mixing is least understood. The function of jet penetration is 
to properly distribute the fuel sources into the combustor by means of jet 
momentum. The turbulence either existing in the airstream or generated by jet 
penetration and combustion is to break up the fuel sources into many small 
lumps stirring up the fuel and air. Finally, mixing is completed by molecular 
diffusion which brings hydrogen and oxygen molecules into intimate contact. 

The molecular diffusion process is slow, therefore, in order to shorten the 
mixing length, it is necessary to generate strong turbulence and small eddies 
at the molecular level to reduce the distance in which the molecules diffuse 
for complete combustion. 

Since turbulent mixing theory is still in an early stage of development, 
it cannot be applied to the solution of this problem at present. Therefore, 
future solution of this problem will have to rely on experiment. 


8.2 FUTURE WORK RECOMMENDATIONS 


The following items are recommended for future work to extend and com- 
pliment present efforts in solving potential problems in supersonic combustion. 

8.2.1 Analytical Approach 

The analytical study may be best summarized by Swithenbank's paper. 
Reference 8-2. In general, the low-wave-number end of the turbulent energy 
spectrum depends on the geometry of the turbulence-producing device, while the 
high-wave number region of the spectrum is almost isotropic and independent 
of the source of eddy formation. The smaller eddies are very important in the 
turbulent mixing process in which two gases are stirred together, to the mole- 
cular level. Since the high-wave-number eddies contain about 80 percent of the 
total turbulent energy, the ability of different systems to achieve mixing 
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will depend largely upon the total turbulent energy In each particular system. 
In general, the ability of a turbulent mixing system to mix two gases will 
depend upon 

(a) Intensity of turbulence 

(b) The scale of turbulence 

(c) Mixing region length 

(d) Lateral distances through which the injected fluid must travel 

(e) Velocity of gas through the mixing system. 

To obtain a semi -quant i tat i ve analysis of the mixing process, we need to 
know the relation between these five various quantities, and the effect of 
the chamber geometry. 

The analytical study should provide guidance in conducting a well defined 
experimental effort. 

8.2.2 Experimental Approach 

The fuel mixing process is strongly influenced by the levels of intensity 
of turbulence. The one-dimensional chemical kinetics program has indicated 
that the separation bubble in the first stage is in the order of 30 to 40 per- 
cent of the duct area, while in the second stage, only 15 percent of flow 
area is blocked. This points to a low level of turbulence, and poor mixing in 
the second' stage. General Electric's Report, Reference l-|, stated: "The 

constant-area section provides maximum pressure disturbances for given fuel 
injection and burning; a necessary situation if high-Mach efficiency is to be 
good . Since the quoted statement was based on experiments, no detailed 
explanation was given in that report. The static pressure distributions from 
this test are reproduced in Figure 8.2-1 for reference. The implication was 
two-fold. First because the "high-Mach" referred to was Mach 9.5, where burner 
inlet static pressure and temperature were approximately 4.8 PSIA and 2500°F 
respectively, the chemical kinetic problem was absent. The calculated ignition 

delay time was at most 50 p.sec. The divergence parameter -i-i- — was 0.0139 
M I A I dX 

cm at the initial point. Second, because of the high static temperature, 
the static pressure rise due to combustion is lower; this results in less flow 
separation, less turbulence and poorer mixing. The latter is considered to 
be the basic problem in future high-Mach number flight conditions. 

In view of the above discussion, the future research program should be 
directed toward the following: 

(a) To correlate turbulence produced by external sources with mixing, 
such as the effects of boundary layer thickness, shock strength, 
size of separation bubble, combustor geometry, pressure gradient and 
combustion. Particular emphasis should be placed on a confined 
geometry with favorable pressure gradients. 
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(b) To develop a turbulence generation method. Swithenbank argued that 
the distance to generate turbulence in a swirling jet is about |/3 

the distance required by a nonswirling Jet. Furthermore, the intensity 
of the turbulence in the developed region is higher in a swirling 
jet, because for a given axial velocity of the jet there is an 
additional tangential velocity representing an increase in the jet 
kinetic energy. 

(c) One of the difficulties in experimental work on turbulent mixing is 
the accurate measurement of the mixing process. The measured quan- 
tities are averages of the instantaneous values because the turbulence is a 
random, fluctuating flow and the flow is three-dimensional. The mixing is not 
as complete as that indicated by the time-averaged measurements taken over a 
period of time using a gas-sampling probe. The sampling probe mixes rich and 
lean regions to indicate an average composition. This may explain the low 
combustor efficiencies in the two-dimensional combustor tests, in which gas- 
sampling measurements indicated uniform mixed flow at the exit. This phenom- 
enon was also reported in Reference 8-3 and requires further exploration. 

8.2.3 Methods to Improve Mixing 

A sudden rearward-facing step in a combustor will promote mixing, how- 
ever, it would result in lower pressure forces on the combustor wall. In 
other words, the rearward-facing step will produce higher chemical efficiency, 
but lower total pressure. The tradeoff between these two should be closely 
examined. 

At high-flight Mach numbers where thermal choking is not likely to occur, 
a converging burner which would enhance flow separation cannot be ruled out as 
a potential solution. 

The effects of surface roughness were clearly demonstrated in the Phase 
1 boilerplate engine. In the paper presented at the 3rd AIAA Propulsion 
Specialist Meeting, it was stated that "The Haste lloy-X-coated engine demon- 
strated quite different characteristics from the pure nickel-coated engine.” 

The high pressure at the fuel injection station indicated spontaneous combus- 
tion, the flow was nearly choked, the ignition temperature was lowered about 
300°F and the slope of the pressure curve was much steeper. This phenomenon 
may be due to a combination of the effects of (l) the higher surface tempera- 
ture, ( 2 ) the increased thickness of Hastelloy X coating, which reduces the 
flow area and increases the static pressure and temperature, ( 3 ) the greater 
area contact with the airstream of the Hastelloy X, which contains 50 percent 
nickel, and ( 4 ) turbulence produced by the rough surface coating. In order to 
further examine these effects, the surface was smoothed from the throat to a 
distance approximately 12 to 14 in. downstream of the throat. The test results 
indicated that there was no significant difference between the rough- and the 
smooth-surfaced engines. The smoothed surface was later verified by laboratory 
check as having an average roughness of 300 microinches. 

Examination of the methods discussed above indicates the surface coating 
to be most adaptable to the AIM engine. 
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9. CONCLUSION 


The results of this study lead to a number of conclusions with regard to 
chemical kinetics in a supersonic combustor. One broad conclusion obtained 
is that chemical kinetics is very sensitive to the divergence of the combustor. 
Basically^ this is attributed to the low stream static temperature which 
"freezes" the species in a manner very similar to that in nozzle flow. 

The detailed chemi ca 1 -ki net i c simulation of the HRE two-dimensional com- 
bustor test results indicated that low combustor efficiencies in the second 
stage were due to mixing. In the extreme case of late but complete mixing, the 
chemical efficiency approached 82 percent which indicated that the second-stage 
flow was not overly kinetics-limited due to the decreasing temperature in the 
divergent duct. 

For the AIM engine at Mach 6 flight conditions, with second-stage fuel 
injected at, stations 2a and 2c, calculations indicated that combustion was 
not chemi ca 1 -ki net i cs limited. Assuming that the fuel was completely mixed 
with air at the injection station, the second-stage flow was choked. 
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RESULTS FROM CHEMICAL KINETIC SIMULATION 
OF RUN 238, TIME - 520.30 SECONDS. 
STOICHIOMETRIC FUEL-AIR RATIO - 0.032568 
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I . Fi rst Stage 
Mixing Model 
Station 12. 


Combustor - Run 238 

2, 50 percent mixing at each injector, 100 percent mixing at 


Station 

CM 

• 

C 

< 

m(lb m/sec) 

m (mixed) 

^2 


c 

V2 

0 

7.3000 

2.450 

1 .7810 

0 

0 

0 


0 

7.5212 

2.518 

1 .7881 

0.003551 

0.50 

0.0071 

0.007078 

0.1220 

7.8112 

2.483 

1 .7881 

0.003551 

0.50 

0.0083 

0.007078 

0.1220 

8.5162 

2.483 

1 .7941 

0.006561 

0.50 

0.0000 

0.013122 

0.2262 

9.0012 

2.434 


0.006830 

0.52 

0.0002 



10.0462 

2.344 


0.009121 

0.70 

0.0041 



10.8462 

2.320 


0.01 1693 

0.89 

0.0443 



11.0912 

2.324 


0.012234 

0.93 

0.0937 



1 1 .4887 

2.390 


0.012839 

0.98 

0.2856 



1 1 .5887 

2.432 


0.012938 

0.99 

0.3668 



1 1 .6887 

2.491 


0.013015 

0.99 

0.5361 



12.0337 

2.907 


0.0I3M9 

1 .00 

0.7912 




12.5337 

13.0212 

13.5737 

14.0287 

14.5112 

15.0062 

15.5562 
16.0562 

16.5562 
17.0500 
17.5500 

17.9950 
18.4950 

18.9950 
19.1950 


3.398 

3.794 

3.970 

4.112 

4.207 

4.305 

4.536 

4.781 

5.067 

5.303 

5.308 

5.253 

5.138 

5.032 

4.992 


I .7941 


0.0131 19 


0.9306 

0.9552 

0.9660 

0.9713 

0.9750 

0.9779 

0.9997 

0.9998 

0.9998 

0.9998 

0.9998 

0.9998 

0.9998 

0.9998 

0.9998 


0.013122 0.2262 
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Figure A.I-I. Temperature vs Axial Distance 
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Figure A. 1-4. Concentration of H^O vs Axial Distance 
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Figure A. 1-9. Concentrations of HO. vs Axial Distance 
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2. Second Stage Combustor - Run 238 

MIxtng Model 2, 50 percent local mixing at each Injector, 100 percent 
mixing at Station 28.62. 

P . m (mixed) 

Station A(ln- ) m(lb m/sec) ”2 


19.1950 

4.990 

1 .7941 

0.0131 19 

19.8550 

5.300 

1 .8068 

0.019491 

20.6000 

5.326 

1 .8068 

0.022315 

20.8100 

5.757 

1.8176 

0.024931 

21 .4150 

6.192 


0.025068 

22.0150 

6.605 


0.025461 

22.5600 

6.943 


0.026041 

23.1600 

7.270 


0.026928 

23.5600 

7.484 


0.027662 

24.1600 

7.818 


0.028977 

24.5600 

8.056 


0.02998 

25.1050 

8.399 


0.03150 

25.5050 

8.631 


0.032516 

26.0750 

8.985 


0.033765 

26.4750 

9.286 


0.034509 

27.0750 

9.755 

1 

0.035420 

27.4750 

10.079 


0.035891 

28.0300 

10.547 

▼ 

0.036363 

28.6200 

1 1 .065 

1.8176 

0.036635 


’m 

’'c 

"’”2 ^ "’”^2 

0 

1 .00 

0.9998 

0.013122 

0.2262 

0.755 

0.6917 

0.025815 

0.4443 

0.755 

0.8612 

0.025815 

0.4443 

0.680 

0.7209 

0.036640 

0.6317 

0.684 

0.8774 



0.695 

0.9150 



0.71 1 

0.9257 



0.735 

0.9284 



0.755 

0.9272 



0.791 

0.9230 



0.819 

0.9191 



0.860 

0.9142 



0.887 

0.9132 



0.922 

0.9135 



0.942 

0.9142 



0.967 

0.9159 



0.980 

0.9175 

1 

1 

0.992 

0.9203 

T 

1 

1 .00 

0.9244 

0.036640 

0.6317 
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Figure A. 2-2. Area Ratio vs Axial Distance 
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Figure A. 2-5. Concentration of OH vs Axial Distance 
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Figure A;2-6. Concentration of 0- vs Axial Distance 















Figure A. 1-7. Concentration of 0 vs Axial Distance 
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Figure A.2-8. Concentration of H vs Axial Distance 
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3. Second Stage Combustor - Run 238 

Mixing Model 3, 0 percent local mixing at injectors, 100 percent mixing at 
Station 29. 


Station 

A(ln.^) 

m(lb m/sec) 

m„ (mixed) 
^2 

’m 


”2 ”^2 

0 

19.1950 

4.990 

1 .7941 

0.013155 

1 .00 

0.9998 

0.0131 15 

0.2261 

19.5437 

4.923 

1 .7941 

0.0131 15 

1 .00 

0.9998 

0.013115 

0.2261 

19.8337 

5.296 

1 .8068 

0.0131 15 

0.508 

0.9872 

0.025815 

0.4451 

20.5687 

5.212 

1 .8068 

0.013115 . 

0.508 

0.9881 

0.025815 

0.4451 

20.8387 

5.549 

I .8176 

0.0131 15 

0.358 

0.9884 

0.036640 

0.6317 

21 .5362 

5.580 


0.013251 

0.362 

0.9770 



22.0812 

5.367 


0.013636 

0.372 

0.9569 



22.4812 

5.725 


0.014050 

0.383 

0.9385 



23.0812 

5.872 


0.014882 

0.406 

0.9084 



23.4812 

5.982 


0.015577 

0.425 

0.8887 



24.0812 

6.172 


0.016828 

0.459 

0.8620 



24.4812 

6.319 


0.017803 

0.486 

0.8468 



25.0312 

6.551 


0.019325 

0.527 

0.8296 



25.6262 

6.795 


0.02121 1 

0.579 

0.8171 



26.0162 

6.985 


0.022581 

0.616 

0.8128 



26.5962 

7.385 


0.024815 

0.677 

0.8103 



26.9962 

7.706 


0.026494 

0.723 

0.8102 



27.5962 

8.263 


0.029221 

0.798 

0.81 19 



27.9962 

8.689 


0.031179 

0.851 

0.8136 



28.3962 

9.157 

1 

0.033249 

0.907 

0.8151 



28.5962 

9.409 


0.34326 

0.937 

0.8157 

1 

1 

28.7962 

9.673 

I 

0.035431 

0.967 

0.8160 

T 

i 

29.0000 

9.955 

1 .8176 

0.036586 

0.999 

0.8162 

0.036640 

0.6317 
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Figure A. 3-11. Concentration of HZ 2 vs Axial Distance 
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4. Second Stage Combustor - Run 238 


Mixing Model 
at Stati on 24 

1 , 50 percent 
.8. 

local mixing 

each 

i njector. 

60 percent 

mi xt ng 

Station 

A(in.^) 

• 

m(lb m/sec) 

m (mixed) 
^2 

TJ 

m 

’Ic 

^2 * 

0 

19.1950 

4.990 

1 .7941 

0.0131 19 

1 .00 

0.9998 

0.013122 

0.2262 

19.5437 

4.923 

1 .7941 

0. 0131 19 

1 .00 

0.9998 

0.013122 

0.2262 

20.0137 

5.273 

1 .8068 

0.015287 

0.59 

0.8697 

0.025815 

0.4451 

20.5687 

5.243 

1 .8068 

0.015287 

0.59 

0.9486 

0.025815 

0.4451 

21 .0287 

5.617 

1 .8176 

0.018020 

0.49 

0.8376 

0.036640 

0.6317 

21 .5612 

5.774 


0.019055 

0.52 

0.8706 



22.1062 

5.987 


0.019951 

0.54 

0.8954 



22.5062 

6.189 


0.020504 

0.56 

0.9070 



23.1062 

6.486 


0.021 167 

0.58 

0.9202 



23.5062 

6.679 


0.021498 

0.59 

0.9272 



24.1062 

6.960 


0.021628 

0.596 

0.9367 



24.5062 

7.144 


0.021938 

0.599 

0.9429 



25.0437 

7.386 


0.021962 

0.599 

0.951 1 



25.6437 

7.595 




0.9570 



26.0212 

7.726 




0.9598 



26.5887 

7.993 




0.9630 



27.1887 

8.295 




0.9658 



27.5887 

8.5t0 

I 



0.9672 



27.9887 

8.739 


1 

1 

0.9684 

1 

1 

28.5887 

9.106 


? 


0.9700 

1 

1 

29.0000 

9.382 

1.8176 

0.021962 

0.599 

0.9709 

0.036640 

0.6317 
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Figure A.4-1. Temperature vs Axial Distance 
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Figure A.4-2. Area Ratio vs Axial Distance 
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Figure A.4-3. Concentration of H„0 vs Axial Distance 
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Figure A. 4-4. Concentration of H- vs Axial Distance 
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Figure A. 4-5. Concentration of OH vs Axial Distance 
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Figure A. 4-7. Concentration of 0 vs Axial Distance 
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Figure A.4-8. Concentration of H vs Axial Distance 
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Figure A. 4-9. Concentration of HO, vs Axial Distance 
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Figure A. 4-10. Concentration of vs Axial Distance 
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Figure A.4-II. Concentration of HZ- vs Axial Distance 
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5. Second Stage Combustor •> Run 238 

Mixing Model I, 50 percent local mixing each injector, 60 percent mixing 
at Station 29. 


Stati on 

A(in.^) 

m(lb m/sec) 

m (mixed) 

”2 



20.8587 

5.607 

1 .8176 

0.017601 

0.48 

0.8304 

21 .5612 

5.764 


0.018322 

0.50 

0.8894 

22.1062 

5.947 


0.018838 

0.51 

0.9122 

22.5062 

6.120 


0.019190 

0.52 

0.921 1 

23.1062 

6.374 


0.019681 

0.54 

0.9290 

23.5062 

6.541 


0.019981 

0.55 

0.9325 

24.1062 

6.794 


0.020392 

0.56 

0.9365 

24.5062 

6.965 


0.020640 

0.56 

0.9386 

25.0437 

7.198 


0.020939 

0.57 

0.9410 

25.6437 

7.408 


0.021229 

0.58 

0.9435 

26.0212 

7.545 


0.021386 

0.583 

0.9451 

26.5887 

7.822 


0.021588 

0.589 

0.9475 

26.9887 

8.028 


0.021705 

0.592 

0.9491 

27.5887 

8.358 


0.021840 

0.596 

0.9516 

28.1887 

8.717 

1 

0.021928 

0.598 

0.9541 

28.5887 

8.973 


0.020960 

0.599 

0.9561 

29.0000 

9.256 

1 .8176 

0.021972 

0.60 

0.9583 



+ m 


Hz. 


0 


0.036640 0.6317 


I 1 
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Figure A. 5-1. Temperature vs Axial Distance 
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Figure A. 5-3. Concentration of H«0 vs Axial Distance 
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Figure A. 5-6. Concentration of 0, vs Axial Distance 
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Figure A. 5-8. Concentration of H vs Axial Distance 
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APPENDIX B 

RESULTS FROM CHEMICAL KINETIC SIMULATION 
OF RUN 242, TIME - 682.55 SECONDS. 
STOICHIOMETRIC FUEL-AIR RATIO - 0.034746 
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I. First Stage Combustor - Run 242 

Mixing Model 2, 50 percent mixing at each injector, iOO percent mixing at 
Station 12. 


Station - ^ m^ (mixed) 

Station A(in. ) m(lb m/sec) "z 



7.3 

3. 126 

1 .536 

7.5125 

3.203 

I .545 

8.0925 

2.974 

1 .545 

8.5925 

2.934 

1 .5524 

9.1925 

2.765 


9.5925 

2.684 


10.1000 

2.796 


10.5000 

3.083 


t 1 .0450 

3.318 


1 1 .4450 

3.190 


1 1 .9525 

3.255 


12.5525 

3.569 


13.0900 

3.885 


13.4900 

4.006 


14.0375 

4.181 


14.6275 

4.259 


15.1975 

4.342 


15.5975 

4.408 


16.1975 

4.512 


16.5975 

4.584 


17.1425 

4.787 


17.5425 

5.071 


18.0800 

5.150 


18.6800 

4.990 

1 

19.0800 

4.893 

T 

19.6000 

4.776 

1 .5524 

9.7925 

2.679 


1 1 .6450 

3.139 


17.7425 

5.234 



0 

0 

0 

0 

0 

0.004415 

0.50 

0 

0.00883 

0. 1654 

0.004415 

0.50 

0 

0.00683 

0. 1654 

0.00821 1 

0.50 

0.0026 

0.01640 

0.3073 

0.008841 

0.54 

0.0353 



0.009797 

0.59 

0.1874 



0.01 1626 

0.71 

0.6259 



0.013387 

0.82 

0.8092 



0.015178 

0.93 

0.8831 



0.015986 

0.97 

0.9155 



0.016396 

1 .0 

0.9436 



0.016400 

1.0 

0.9586 




0.9770 
0.9673 
0.9704 
0.9729 
0.9749 
0.9761 
0.9776 
0.9784 
0.9794 
0.9798 
0.9801 
0.9809 
0.9815 

0.016400 1.0 0.9824 0.01640 0.3073 
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Figure B.I-IO. Concentration of H 2 O 2 vs Axial Distance 
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2* Second Stage Combustor - Run 242 

Mixing Model 2, 50 percent mixing at each Injector, 100 percent mixing at 
Station 28.75. 


Stati on 

A(in.^) 

m(lb m/sec) 

m (mfxed) 

”2 

n 

m 


^2 ^ V 2 

0 

19.6000 

4.762 

1 .55240 

0.016400 

1 .0 

0.9824 

0.01640 

0.3073 

19.8262 

5.443 

1 .57318 

0.026800 

0.27 

0.6738 

0.03718 

0.6966 

20.5112 

5.648 

1 .57318 

0.026800 

0.27 

0.8854 

0.03718 

0.6966 

20.9012 

6.508 

1 .59100 

0.035730 

0.64 

0.821 1 

0.05500 

1 . 0305 

21 .5012 

6.915 


0.036007 

0.65 

0.8831 



22.0237 

7.093 


0.036584 

0.67 

0.9043 



22.6062 

7.509 


0.037597 

0.68 

0.9129 



23.0737 

7.815 


0.038692 

0.70 

0.9118 



23.6737 

7.904 


0.040465 

0.74 

0.9098 



24.2t 12 

8.224 


0.042404 

0.77 

0.8999 



24.BI 12 

9.269 


0.044961 

0.82 

0.8846 



25.5487 

9.426 


0.048187 

0.88 

0.8657 



26.2862 

9.503 


0.050797 

0.92 

0.8586 



27.0287 

10.169 


0.052794 

0.96 

0.8537 


1 

27.8112 

9.538 


0.054213 

0.99 

0.8539 

1 


28.5487 

10.065 


0.054908 

0.999 

0.8607 

? 


28.7487 

10.334 

1 .59100 

0.055000 

1 .0 

0.8624 

0.05500 

1 .0305 

19.6912 

5.043 







19.9762 

5.451 







25.0062 

9.725 







25.9487 

9.278 







28.0062 

9.448 
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MIXING MODEL 2 



Figure B.2-2. Area Ratio vs Axial Distance 
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Figure B.2-II. Concentration of HZ, vs Axial Distance 
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Figure B.2-5. Concentration of OH vs Axial Distance 
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Figure B.2-6. Concentration of 0 vs Axial Distance 
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Figure B. 2-7. Concent rat i on of 0^ vs Axial 0 i stance 
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Figure B.2-10. Concentration of H,0, vs Axial Distance 
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3. Second Stage Combustor - Run 242 

Mixing Model I, 50 percent local mixing at each injector, 60 percent 
mixing at Station 24.9. 


2 . ">0 (mixed) 
Station A(in. ) m(lb m/sec) "2 

19.6000 4.762 
19.8825 5.428 
20.4825 5.390 
20.9550 6.048 
21.5150 6.218 
22.0600 6.399 
22.6600 6.852 
23.0812 7.137 
23.6462 7.223 
24.0237 7.290 
24.5912 8.072 
24.9912 8.861 
25.7287 8.426 
26.0662 8.316 
26.6662 8.715 
27.0287 8.943 
27.6112 8.357 
28.0062 8.073 
28.5487 8.586 
28.7487 8.816 
29.0000 9.144 


’^m 

’Ic 

”2 ”^2 

0 

1 .0 

0.9824 

0.01640 

0.3073 

0.55 

0.8382 

0.03718 

0.6966 

0.55 

0.8382 

0.03718 

0.6966 

0.46 

0.8567 

0.05500 

1 .0305 

0.49 

0.8869 



0.52 

0.9002 



0.55 

0.9086 



0.57 

0.9120 



0.58 

0.9181 



0.59 

0.9234 



0.596 

0.9302 



0.597 

0.9344 




0.9400 




0.9431 




0.9476 




0.9495 




0.9526 



t 

0.9554 

1 



0.9587 


( 

\ 

0.9594 

» 

i 

0.597 

0.9602 

0.05500 

1 .0305 


1.5524 0.016400 

I .5732 0.020382 

I .5732 0.020382 

1.5910 0.025140 

0.027167 
0.028842 
0.030347 
0.031 191 
0.032048 
0.032445 
0.032777 
0.032824 


1 i 

1.5910 0.032824 
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Figure B.3-3. Concentration of H^O vs Axial Distance 
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Figure B.3-4. Concentration of H_ vs Axial Distance 
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Figure B.3-5. Concentration of OH vs Axial Distance 
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Figure B.3-7. Concentration of 0 vs Axial Distance 
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Figure B.3-8. Concentration of H vs Axiai Distance 
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Figure B.3-10. Concentration of H»0. vs Axial Distance 
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Figure B.3-II. Concentration of HZ, vs Axial Distance 
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4. Second Stage Combustor - Run 242 

Mixing Model I, 50 percent local mixing at each injector, 60 percent 
mixing at Station 29. 



9 


m (mixed) 

^2 





Stati on 

A(in. ) 

m(lb m/sec) 

'^m 

^c 

”2 »"2 


19.6000 

4.762 

1 .5524 

0.016400 

1.0 

0.9824 

0.01640 

0.3073 

20.4825 

5.390 

1 .5732 

0.020382 

0.55 

0.9232 

0.03718 

0.6966 

20.9050 

6.036 

1 .5910 

0.024745 

0.45 

0.8547 

0.05500 

1 .0305 

21 .2162 

6.114 

1 .5910 

0.025355 

0.46 

0.8877 



21 .6162 

6.206 


0.026104 

0.47 

0.9024 



22.3537 

6.507 


0.027381 

0.50 

0.9181 



22.9287 

6.900 


0.028284 

0.51 

0.9223 



23.6312 

7.020 


0.029276 

0.53 

0.9256 



24.3762 

7.519 


0.030195 

0.55 

0.9292 



24.9762 

8.588 


0.030836 

0.56 

0.9279 



25.1137 

8.587 


0.030970 

0.56 

0.9271 



25.3137 

8.457 


0.031 157 

0.57 

0.9265 



25.5137 

8.339 


0.031334 

0.57 

0.9266 



25.7137 

8.231 


0.031501 

0.57 

0.9271 



26.4512 

8.394 


0.032032 

0.58 

0.9304 



27 . 1 962 

8.599 


0.032432 

0.59 

0.9326 



27.9962 

7.964 

1 

0.032710 

0.59 

0.9388 



28.5337 

8.474 


0.032807 

0.596 

0.9437 

1 

1 

28.7337 

8.709 


0.032825 

0.597 

0.9451 

f 

? 

29.0000 

9.060 

1 .5910 

0.032834 

0.597 

0.9467 

0.05500 

1 . 0305 


19.7975 5.443 
22.0162 6.282 
24.0162 7.071 
24.1762 7.257 
24.5762 7.821 
24.7762 8.172 
26.0512 8.122 
26.2512 8.255 
27.0212 8.790 
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Figure B.4-2. Area Ratio vs Axial Distance 
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Figure B.4-3. Concentration of H,0 vs Axial Distance 
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Figure B.4-4. Concentration of H, vs Axial Distance 
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Figure B.4-5. Concentration of OH vs Axial Distance 
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Figure B.4-6. Concentration of 0^ vs Axial Distance 
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Figure B.4-7. Concentration of 0 vs Axial Distance 
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Figure B.4-8. Concentration of H vs Axial Distance 
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Figure B.4-9. Concentration of HO, vs Axial Distance 
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Figure B.4-10. Concentration of H-O- vs Axial Distance 



















